NAVAL  POSTGRADUATE  SCHOOL 
Monterey ,  California 


THESIS 


A  DAY  IN  Till:  LIFE  OF  A  WARM  FROM 

by 

Susan  A.  Davies 

December  1989 


Thesis  Advisor 


Wendell  A.  Nuss 


Approved  for  public  release;  distribution  is  unlimited. 


REPORT  DOCUMENTATION  PAGE 

la  Report  Sccu:  t\  Cias?i':ca::on  Unclassified 

ib  Restrictive  Markings 

2aSecuni:  Clas    ficati  m  \uthorin 

3  Distribution  Availability  of  Report 

Approved  for  public  release:  distribution  is  unlimited. 

IPe      rmiru n  Repori  Number(s) 

5  Monitoring  Or-c.n:;';-/:  ':    K.;;-.-'-  \umbensi 

Naval  Postgraduate  School 

6b  Office  Symbol 
(if  applicable)  i5 

7a  Name  of  Monitoring  Organization 
Naval  Postgraduate  School 

Monterey.  CA  93943-5000 

7b  Address  i  city,  state,  and  ZIP  code) 
Monterey,  CA  93943-5000 

Si  Name  of  Furc:::;  Sponsoring  Organization      |  Sb  Office  Symbol 

9  Procurement  Instrument  Identification  Number 

jic  \ddre  >      .;..    rra«    and  ZIP  code) 

10  Source  of  Fundina  Numbers 

Program  Element  No  |  Project  No   |  Task  No   ]  Work  I  nit  Accession  No 

|1  Title    include  security    la    Iflcathr. )  A  DAY  IN  THE  LIFE  OF 

A  WARM  FRONT 

2]     sonaiAui          Susan  A.  Davies 

fea  Typed  Rep  ■:: 
Master's  Thesis 

1       .  Covered 
From                       To 

14  Date  of  Report  i  year,  month,  day)             1  15  Page  Count 
December  1989                                   |S1 

Sup;  ji-entar;  Nota:       rhe  views  expressed  in  this  tl     i    an 
ption  of  the  Department  of  Defense  or  the  U.S   G  »   n  mei 

those  of  the  author  and  do  not  reflect  the  official  policy  or  po- 

IT  C<  sat   Codes 

Warm  front  secondary  cir 

everse  if  r*                                                 numl  r',1 
culation,  influence  of  surface  fluxes 

G       -             S      :     .:■ 

iributk 


aluate  the  effect  on  t 


Three  numerical  model  simulation 
nhancing  the  comp  tnents  forcing  the 
fcuation  is  mine  the  cc 

pibutions  are  varied  in  each  case,  in  oi 

Results  confirm  previ   j^  studies  thai  geostrophic  deformation  forces  stro; 

^ontogenesis  is  weaker  and  forced  primarih  by  latent  heat  release.    Frontogi 

tional  forcing  and  diabatic  heating,  v.  huh  depend  upon  the  surface  and  bound 

five  percent  ol  total  forcing  of  the  warm  front,  it  apparen 

prcct  effeel  re  dominant  processes    The  intensity  of  frictional  fc 


are  compared  to  examine  the  role  of  boundary  layer  stratification  in 


c  circulation  of  an  idealized  warm  front.    The  Sawyer-Eli 
ol  frictional  forcing  and  diabatic  heatins 


m  Iroi 
ompor 


n  diagnostii 
as  well  as  confluent  and  shear 
ce  heat  and  moisture  flux  dis- 
-  the  cyclogenesis. 
it  the  surface,  and  at  mid-levels 
is  modified  by  small-scale  fric- 
■es.  Although  frictional  forcing 
ntogenesis,  partialis  due  to  in- 
dependent on  the  surface  heat 
iux  distribution  and  track  of  the  cyclone  relative  to  the  sea-surface  temperature  gradient.  In  the  absence  of  surface  fluxes. 
able  Surface  forcing  due  to  diabatic  heating  is  frontolytical,  and  reduces  total  forcing  at  low  levels 
by  about  twenty  percenl  The  magnitude  of  the  frontolytical  forcing  by  diabatic  heating  at  low  levels  is  only  partially  de- 
pendent on  the  surface  heat  and  moisture  flux  distribution.  Even  in  the  absence  of  surface  heat  and  moisture  fluxes,  the 
m  ntolytical  forcing  persisted,  suggesting  problems  with  the  cumulus  parameterization  or  errors  in  the  computation  of  this 
term.    These  results  remain  to  be  verified  32ainst  an  actual  case  studv  of  a  warm  front. 


'I  he  sur 
1  as  well 


er  proce 
lances  h 
s  strong 


I   Abstract 

Z     arm  a 


21  Abstract  Security  Clas 
I  'nclassified 


Wendell  A.  Nuss 


[408)  646-3275 


22c  Office  Symbol 

63NU 


ID  i  ORM  1473,84  MAR 


y  classification  of  this  page 


Unclassified 


T245413 


Approved  for  public  release:  distribution  is  unlimited. 

A  Day  in  the  Life  of  a  Warm  Front 

by 

Susan  A.  Davies 

/ 

Lieutenant  Commander.  United  States  Navy 

B.A..  University  of  California.  San  Diego.  1980 

Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  METEOROLOGY  AND  PHYSICAL 
OCEANOGRAPHY 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 

December  19S9 


ABSTRACT 

Three  numerical  model  simulations  of  cyclogenesis  are  compared  to  examine  the  role 
of  boundary  layer  stratification  in  enhancing  the  components  forcing  the  ageostrophic 
circulation  of  an  idealized  warm  front.  The  Sawyer-Eliassen  diagnostic  equation  is  ap- 
plied to  examine  the  contributions  of  frictional  forcing  and  diabatic  heating,  as  well  as 
confluent  and  shear  geostrophic  deformation,  to  forcing  of  the  secondary  circulation  of 
a  warm  front.  The  surface  heat  and  moisture  flux  distributions  are  varied  in  each  case, 
in  order  to  evaluate  the  effect  on  each  component  as  well  as  the  cyclogenesis. 

Results  confirm  previous  studies  that  geostrophic  deformation  forces  strong 
frontogenesis  at  the  surface,  and  at  mid-levels  frontogenesis  is  weaker  and  forced  pri- 
marih  by  latent  heat  release,  rrontogenetkal  forcing  is  modified  by  small-scale  fric- 
tional forcing  and  diabatic  heating,  which  depend  upon  the  surface  and  boundary  layer 
processes.  Although  frictional  forcing  comprises  less  than  five  percent  of  total  forcing 
of  the  warm  front,  it  apparently  enhances  frontogenesis.  partially  due  to  indirect  effects 
on  other  more  dominant  processes.  The  intensity  of  frictional  forcing  is  strong!}"  de- 
pendent on  the  surface  heat  flux  distribution  and  track  of  the  cyclone  relative  to  the 
sea-surface  temperature  gradient.  In  the  absence  of  surface  fluxes,  frictional  forcing  is 
_  Surface  forcing  due  to  diabatic  heating  is  frontolytical.  and  reduces  total 
'  low  levels  by  about  twenty  percent.  The  magnitude  of  the  frontolytical  forcing 
by  diabatic  heating  at  low  levels  is  only  partially  dependent  on  the  surface  heat  and 
moisture  flux  distribution.  Even  in  the  absence  of  surface  heat  and  moisture  fluxes,  the 
frontolytical  forcing  persisted,  suggesting  problems  with  the  cumulus  parameterization 
or  errors  in  the  computation  of  this  term.  These  results  remain  to  be  verified  against  an 
actual  case  studv  of  a  warm  front. 


C.i 

TABLE  OF  CONTENTS 

I.  INTRODUCTION    1 

A.  CYCLOGENESIS  AND  FRONTOGENESIS    1 

1.  Sub-Synoptic  Scale  Processes     1 

2.  Planetary 'Boundary  Layer  Processes     2 

B.  MARINE  CYCLONES     3 

C.  OBJECTIVE     5 

II.  BACKGROUND     6 

A.  AGEOSTROPHIC  CIRCULATION  AROUND  A  WARM  FRONT     6 

B.  SAWYER-ELIASSEN  EQUATION     11 

C.  DESCRIPTION  OF  TERMS     13 

1.  Gcostrophic  Confluence     13 

2.  Gcostrophic  Shear     14 

3.  Friction     15 

-4.     Diabatic  1  leating     16 

III.  DESCRIPI  ION  OE  DATA  SET    IS 

A.  NUMERICAL  MODEI 18 

B.  INITIAL  CONDITIONS     VI 

C.  DEVELOPMENT  OE  "THE  CYCLONE7     22 

D.  1  OTAL  FORCING  AT  THE  SURFACE    25 

IV.  RESULTS    30 

A.  OVERVIEW     30 

B.  TOTAL  FORCING  OF  THE  WARM  FRONT   30 

1.  Surface    30 

2.  Cross  Sections      31 

C.  I  OTAL  FORCING  AT  24  HOURS    36 

D.  CONTRIBUTIONS  BY  EACH  COMPONENT  OF  FORCING     37 

1 .  Oven  iew     37 

2.  Case  1     37 


a.  Surface     38 

b.  Cross-Sections    41 

3.  Case  2     43 

a.  Surface     43 

b.  Cross-Sections    46 

4.  Case  3     49 

a.  Surface     49 

b.  Cross-Sections    52 

E.  TIME  SERIES  OF  FORCING    54 

1.  Boundary  Layer     55 

2.  Mid  and  Upper  Levels     57 

F.  COMPARISON  OF  CASES    61 

1 .  Boundary  Layer     61 

2.  Mid  and  Upper  Levels     61 

V.  DISCUSSION    62 

A.  OVERVIEW     62 

B.  BOUNDARY  LAYER     62 

1.  Geostrophic  Deformational  Forcing     62 

2.  Frictional  Forcing    63 

3.  Diabatic  1 1  eating     64 

C.  MID  AND  UPPER  LEVELS     64 

1.  Geostrophic  Deformation    64 

2.  Diabatic  1  leating     65 

VI.  CONCLUSIONS    66 

VIE     RECOMMENDATIONS    67 

LIST  OF  REFERENCES    68 

INI  I IAL  DISTRIBUTION  LIST     71 


LIST  OF  FIGURES 

1.  Streamlines  around  a  positive  point  source 7 

2.  Idealized  cross  section  through  a  warm  front 8 

3.  Secondary  circulation  around  a  warm  front 9 

4.  Agcostrophic  forcing  centers  around  a  warm  front 10 

5.  Transverse  ageostrophic  circulations  for  an  idealized  front 14 

6.  Centers  of  frictional  forcing 16 

7.  Centers  of  diabatic  forcing 17 

8.  Initial  conditions 20 

9.  Case  1  and  Case  2  SST  distribution 21 

10.  As  in  Fig.  9.  except  for  Case  3 22 

11.  SLP  (mb)  versus  time  for  each  case 23 

12.  Total  forcing  at  the  surface  at  24  h.  Case  1 27 

13.  As  in  Fig.  12.  except  for  Case  2 2S 

14.  As  in  Fig.  12.  except  for  Case  3 29 

15.  F  volution  of  total  forcing  of  the  surface  warm  front 31 

16.  24-h  cross  section  of  total  ageostrophic  forcing.  Case  1 32 

]".    As  in  Fig.  16.  except  for  Case  2 33 

18.  As  in  Fig.  16.  except  for  Case  3 34 

19.  Evolution  of  total  forcing  in  the  PDF 35 

2' i.     Evolution  of  total  forcing  at  mid  levels 36 

21.  24-h  heat  flux  distributions  for  Case  1 38 

22.  24-h  surface  forcing  due  to  friction.  Case  1 39 

23.  24-h  surface  forcing  due  to  diabatic  heating.  Case  1 40 

24.  24-h  cros^-section  of  frictional  forcing.  Case  1 41 

25.  As  in  Fig.  24.  except  for  diabatic  heating 42 

26.  24-h  confluent  geostrophic  forcing  at  the  surface.  Case  2 44 

2".    As  in  Fig.  26.  except  for  forcing  due  to  shear 45 

28.  As  in  Fig.  26,  except  for  diabatic  heating 46 

29.  24-h  cross-section  of  confluent  geostrophic  forcing.  Case  2 47 

3').    As  in  Fig.  29,  except  for  forcing  due  to  shear 48 

31.    As  in  Fig.  29.  except  for  diabatic  heating 49 


Fie.  32 


1  '& 

Fig 

33 

Fig 

34 

Fig 

35 

Fig 

36 

Fig 

37 

Fig 

3S 

Fig 

39 

Fie 

40 

Fig 

41 

Fia 

42 

As  in  Fig.  21.  except  for  Case  3 50 

24-h  surface  forcing  due  to  friction.  Case  3 51 

24-h  surface  forcing  due  to  diabatic  heating.  Case  3 52 

24-h  cross-section  of  frictional  forcing,  Case  3 53 

As  in  Fig.  35.  except  for  diabatic  heating 54 

Time  series  of  forcing  in  the  PBL 56 

As  in  Fig.  37,  except  for  Case  3 57 

As  in  Fig.  37.  except  for  600  mb 58 

As  in  Fig.  38.  except  for  600  mb 59 

As  in  Fig.  37.  except  for  300  mb 60 

As  in  Fig.  38.  except  for  300  mb 61 


ACKNOWLEDGEMENT 

My  thanks  to  Professor  Wendell  Nuss  for  his  patience  with  an  aviator. 


I.     INTRODUCTION 

A.     CVCLOGENESIS  AND  FRONTOGENESIS 

Mid-latitude  cyclones  are  transient  waves  in  the  global  zonal  flow  that  perform  an 
important  role  in  transporting  heat  and  momentum  in  the  atmosphere.  Historically, 
cyclone  development  was  explained  as  a  wave  instability  on  the  polar  front,  with  the 
existence  of  the  front  considered  a  necessary  precursor  for  cyclogenesis.  This  theory. 
known  as  the  polar  front  theory  of  cyclogenesis.  significantly  influenced  Petterssen 
(1956)  who  related  cyclogenesis  to  the  interaction  of  vorticity  advection  aloft  with  sur- 
face frontal  features.  The  modern  view,  based  on  baroclinic  instability  theory,  considers 
fronts  to  be  a  consequence  of  the  development  of  a  cyclone  and  a  source  of  baroclinic 
energy  for  cyclogenesis.  The  baroclinicity  of  the  atmosphere  and  the  synoptic-scale  flow 
combine  to  lead  to  cyclogenesis.  the  development  of  sub-synoptic  scale  fronts,  and  the 
tion  of  jet  streaks. 
The  basic  premise  of  baroclinic  instability  theory  is  that  the  circulation  of  a  cyclone 
re  ilts  from  an}  small  perturbation  forced  by  advection  of  positive  (cyclonic)  vorticity 
aloft  with  the  advance  of  the  upper-level  wave  or  by  warm  advection  at  low  levels  over 
a  broad  region.  Ge  istrophic  deformation  associated  with  cyclogenesis  acts  to  enhance 
the  temperature  gradient  in  particular  regions  of  a  cyclone.  This  process  is  known  as 
mesis.  Vertical  and  divergent  ageostrophic  motions  arise  in  response  to 
frontogenesjs.  arid  are  the  mechanism  for  the  adjustment  of  vorticity  and  thermal  prop- 
erties to  return  the  system  to  thermal  wind  balance.  The  ageostrophic  circulations  as- 
sociated with  frontogenesis  are  small-scale  manifestations  of  the  basic  synoptic-scale 
forcing.  Although  these  synoptic-scale  processes  are  thought  to  be  the  primary  proc- 
esses for  cyclogenesis,  frontogenesis  is  a  secondary  process  that  may  feed  back  on  the 
synoptic-scale  forcing  to  produce  stronger  cyclogenesis  than  that  given  by  just  the 
synoptic-scale  forcing.  The  interaction  between  synoptic-scale  vorticity  and  temperature 
advection  and  mesoscale  frontogenesis  in  cyclones  is  not  yet  completely  understood. 

1.     Sub-Synoptic  Scale  Processes 

Recent  c\clogenesis  studies  suggest  important  coupling  between  cyclogenesis 
and  mesoscale  processes  associated  with  frontogenesis  and  jet  streaks.  L'ccellini  et  al. 
(1985,  1987)  and  Keyser  and  Shapiro  (19So)  have  shown  that  small  scale  upper-level 
baroclinic  processes  associated  with  upper-level  jet  streaks  enhance  the  vertical  circu- 


lation  and  cyclogenesis.  Strong  jet  streak  divergence  induces  a  secondary  circulation 
that  extends  through  the  depth  of  the  troposphere.  This  enhances  the  development  of 
a  low-level  jet  and  cyclogenesis  at  the  surface.  This  mesoscale  process  is  considered  to 
be  an  important  mechanism  for  some  explosive  cyclones. 

Another  important  mechanism  in  some  explosive  cyclones  is  the  release  of  latent 
heat  in  the  mid-troposphere.  The  release  of  latent  heat  also  influences  th«:  If  ontogenesis 
process.  Williams  et -al.  (1981).  using  a  numerical  model,  and  Thorpe  \19S4).  using  an 
analytical  quasi-geostrophic  model,  showed  that  latent  heat  release  on  the  warm  side  of 
the  front  is  frontogenetical  aloft.  The  heating  directly  increases  the  temperature  gradient 
at  levels  of  large  heating,  and  the  induced  divergent  circulation  leads  to  convergence 
below  the  heating  which  increases  the  temperature  gradient  in  that  region.  The  precipi- 
tation maximum  occur';  where  the  vertical  motion  is  a  maximum.  Evaporative  cooling 
occurs  when  precipitation  falls  through  dry  air  below  on  the  cold  side  of  fronts,  thus 
enhancing  the  cross-front  temperature  gradient  ahead  of  the  front. 
2.     Planetary  Boundary  Layer  Processes 

Although  baroclinic  synoptic  and  mesoscale  processes  mentioned  above  arc 
dominant,  the  importance  of  surface  and  boundary  layer  processes  also  has  been  dem- 
onstrated in  recent  studies.  The  boundary  layer  processes  contribute  to  cyclogenesis  and 
frontogenesis  through  thermodynamic  effects  due  to  the  heat  and  moisture  lluxes.  and 
dynamic  effects  due  to  friction.  Friction  further  disrupts  geostrophic  balance  and  in- 
creases the  ageostrophic  component  of  the  wind  through  frictional  wind  turning.  This 
results  in  an  additional  component  of  convergence  and  shear  due  to  friction.  Keyser  and 
Anthes  ( 19S2)  showed  that  convergence  due  to  surface  friction  acts  frontogenetically  at 
the  surface.  Behind  the  frontal  zone  the  vertical  shear  of  the  ageostrophic  wind,  which 
is  frictionally  induced,  acts  to  diminish  the  static  stability  in  the  planetary  boundary 
layer  (PBL)  and  enhance  the  static  stability  at  the  top  of  the  boundary  layer.  Ahead  of 
the  front,  differential  temperature  advection  by  the  ageostrophic  component  of  the  wind 
acts  to  increase  the  static  stability  in  the  boundary  layer  and  decrease  the  static  stability 
at  the  top  of  the  boundary  layer.  The  result  is  a  compressed  thermal  gradient  in  the 
PBL.  winch  is  frontogenetical.  Keyser  and  Anthes  <  1 CJ S 2 )  also  showed  that  frictional 
convergence  drives  a  vertical  velocity  jet  above  the  surface  front.  This  vertical  jet  is  a 
component  of  the  secondary  circulation. 

Danard  and  Ellenton  (1980)  evaluated  the  surface  heating  in  several  east  coast 
L*.  S.  cyclones,  and  concluded  that  the  surface  heating  acted  to  decrease  the  low-level 
baroclinicin   bv  warming  the  cold  advection  regions  and  vice  versa.     '1  his  finding  is 


consistent  with  previous  studies  of  the  effects  of  surface  heating  in  baroclinic  instability 
models,  and  indicates  a  frontolytical  role  of  the  surface  heat  flux.  However,  in  studies 
of  the  Presidents'  Day  storm  of  1979.  Bosart  (1981)  and  Bosart  and  Lin  (1984)  indicated 
that  the  low-level  contributions  by  the  surface  heat  and  moisture  fluxes  were  important 
to  cyclogenesis.  1  hey  found  that  positive  heat  fluxes  to  the  northeast  of  the  surface  low 
potentially  contributed  to  rapid  development  through  distribution  of  the  cyclogenetic 
environment.  Reed  and  Albright  (1986)  found  similar  results  in  their  study  of  an  eastern 
North  Pacific  cyclone. 

A  favorable  surface  heat  and  moisture  flux  distribution  can  also  enhance 
frontogenesis  and  cyclogenesis  through  the  associated  latent  heat  release  above  the 
boundary  layer.  Chen  et  al.  (1983)  found  that  mid-tropospheric  latent  heating  decreased 
when  the  surface  heat  and  moisture  fluxes  were  removed.  Atlas  (19S7)  found  that  the 
vertical  transport  of  air  that  is  heated  and  moistened  by  the  surface  heat  and  moisture 
fluxes  can  result  in  significant  mid-tropospheric  heating  through  the  release  of  latent 
he  it. 

B        MARINE  CYCLONES 

Compared  with  continental  and  coastal  cyclones,  marine  cyclones  have  a  more 
surface  heat  and  moisture  flux  distribution  contributing  to  frontogenesis  and 
cyclogenesis.  due  to  the  influence  of  the  sea-surface  temperature  (SST)  distribution. 
Businger  and  Shaw  (1984)  suggested  that  for  a  simple  uniform  flow  across  an  SST  gra- 
dient. PHI.  convergence  will  result  from  stratification  effects  on  the  frictional  wind 
is  the  wind  flows  across  an  SST  gradient.  As  a  cyclone  moves  across  a  fixed 
SST  distribution,  changes  in  the  frictional  convergence  and  the  transport  of  boundary 
layer  air  to  the  mid-troposphere  are  likely  to  occur.  In  their  study  of  surface  fluxes  and 
divergence  in  Storm  Transfer  and  Response  Experiment  (STREX)  storms.  Fleagle  and 
Nuss  (19S5)  showed  that  the  distributions  of  the  frictional  convergence  and  vertical  ve- 
locity at  the  top  of  the  PBL  are  primarily  determined  by  the  horizontal  distribution  of 
the  surface  wind  stress.  In  turn,  the  wind  stress  depends  not  only  on  the  synoptic-scale 
forcing  at  the  top  of  the  PBL.  but  also  on  the  PBL  stratification,  the  surface  heat  flux 
distribution,  and  thermal  wind  turning  in  the  PBL.  Nuss  (1989)  showed  that  stable 
stratification  behind  the  warm  front,  and  unstable  stratification  ahead  of  the  warm  front 
in  a  region  of  upward  heat  fluxes,  apparently  resulted  in  enhanced  frictional  convergence 
ahead  of  the  front.   Danard  (1986)  found  that  increasing  the  SST  by  five  degrees  resulted 


in  greater  boundary  layer  convergence  and  enhanced  convective  heating  in  the  model 
cyclogenesis. 

Nuss  and  Anthes  (1987)  studied  the  effect  of  surface  fluxes  and  boundary  layer  in- 
fluences on  an  idealized  marine  cyclone.  They  showed  that  although  baroclinic  insta- 
bility, which  depends  on  the  vertical  wind  shear  and  static  stability,  as  shown  originally 
by  Charney  (1947)  and  Eady  (1949).  is  sufficient  for  rapid  cyclogenesis.  relatively  small 
changes  in  the  low-level  static  stability,  meridional  temperature  gradient,  and  surface 
heat  and  moisture  fluxes  significantly  modified  the  development  of  the  model  cyclone. 
These  low-level  changes  resulted  entirely  from  varying  the  SST  distribution.  They  found 
that  the  deepening  rate  of  the  cyclone  increased  by  25%  for  complete  removal  of  the 
surface  heat  and  moisture  fluxes,  and  increased  by  15%  for  a  modified  SST  distribution 
in  phase  with  the  low-level  temperature  pattern. 

These  results  were  further  analyzed  by  Nuss  (1989),  in  order  to  clarify  the  role  of 
surface  fluxes  in  the  overall  development  of  the  cyclone.  He  found  that  the  effect  of 
surface  heat  and  moisture  fluxes  on  the  structure  and  dynamics  of  oceanic  cyclones  can 
be  substantial  and  that  strictly  boundary  layer  processes  act  to  enhance  the  vertical  cir- 
culation and  cyclogenesis  through  a  dynamic  feedback  process.  Specifically,  he  found 
that  the  effect  of  surface  heat  and  moisture  fluxes  depends  critically  on  their  distribution 
relative  to  the  cyclone.  For  a  zonal  SST  distribution,  patterns  of  surface  heating 
counteract  PBL  temperature  advection  and  oppose  cyclogenesis.  as  found  in  previous 
studies,  including  Danard  and  Ellenton  (1980).  For  a  sinusoidal  SST  distribution  es- 
sentially in  phase  with  the  low-level  temperature  distribution,  the  resultant  pattern  of 
surface  heating  yielded  unstable  stratification  to  the  northeast  and  stable  stratification 
south  of  the  warm  front.  This  acted  to  enhance  cyclogenesis  through  enhanced  frictional 
convergence  as  well  as  increased  heat  and  moisture  transport  to  the  surface  cyclone  from 
the  northeast. 

The  relationship  of  enhanced  frictional  convergence  to  other  frontogencsis  processes 
was  not  examined  and  is  the  focus  of  this  study.  Although  the  surface  heating  pattern 
may  be  a  damping  effect,  associated  mid-tropospheric  latent  heat  release  increases  the 
temperature  gradient  and  cyclogenesis.  which  indirectly  enhances  the  PBL  frontogencsis 
processes.  '1  his  indirect  relationship  between  surface  diabatic  processes  and  the  intensity 
of  the  geostrophic  frontogencsis  and  cyclogenesis  is  not  well  understood.  Uccellini  et 
al.  (19S7)  have  demonstrated  the  importance  of  the  non-linear  coupling  between 
boundary  layer  processes,  latent  heat  release,  and  jet  stream  forcing  along  a  coastal  front 
foi  the  Presidents'  I)a\  storm.    For  open  ocean  cyclones,  the  warm  front  region  is  likely 


to  be  the  focal  point  of  the  interaction  between  surface  fluxes,  diabatic  heating,  and  the 
jet.  as  suggested  by  Nuss  (19S9).  and  will  be  the  subject  of  this  study.  Studies  by 
Emanuel  ( 1985,  1()SS)  and  Sanders  (1986)  also  suggest  that  latent  heat  processes  associ- 
ated with  the  warm  front  may  contribute  to  enhanced  cyclogenesis  over  the  ocean. 

C.     OBJECTIVE 

The  purpose  of  this  research  is  to  quantify  the  contributions  of  various  forcing 
mechanisms  to  the  secondary  circulation  around  a  warm  front,  and  assess  their  impact 
on  the  development  of  a  cyclone.  The  influence  of  the  surface  heat  and  moisture  fluxes 
on  the  frictional  forcing  and  diabatic  heating,  and  resultant  secondary  circulation,  is  of 
particular  interest  since  it  is  not  as  clearly  understood  as  the  larger-scale  gcostrophic 
deformation  processes.  In  addition,  the  feedback  of  these  diabatic  processes  on  the 
_-  trophic  deformation  and  warm  frontal  intensification  will  be  examined.  The  specific 
objectives  are: 

1.  Diagnose  the  geostrophic  versus  diabatic  processes  in  warm  frontogenesis 

to  distinguish  frictional  convergence  from  pur    geostrophic  deformation; 

2.  Describe  the  feedback  mechanism  of  surface  heat  fluxes  on  geostrophic  and 

diabatic  forcing  by  quantifying  their  Irontogenetic  effects  on  the  friction 
and  latent  heat  release  and  their  frontolytical  effects  on  the  PBL  thermal 
gradient: 

3.  Describe  the  time  evolution  of  the  forcing; 

•4.    Obtain  an  understanding  of  the  interaction  of  the  SST  with  fronts. 

The  Sawyer-Eliassen  diagnostic  equation  will  be  used  to  diagnose  the  forcing  of  the 
vertical  circulation  around  the  warm  front  in  three  numerical  simulations  of 
cyclogenesis.  The  Sawyer-Eliassen  equation  is  reviewed  in  Section  II.  followed  by  a  brief 
description  o[  the  model  experiments  in  Section  III.  Section  IV  compares  the  results 
of  the  three  simulations.  A  discussion  of  these  results  is  given  in  Section  V,  followed  by 
conclusions  and  recommendations  for  future  research  in  Sections  VI  and  VIE 


II.     BACKGROUND 

Sawyer  (1956)  showed  how  the  secondary  circulation  around  fronts  can  be  calcu- 
lated on  the  basis  of  quasi-geostrophic  theory,  and  Eliassen  (1959,  1962)  extended  his 
work  to  include  friction  and  latent  heating  effects,  as  well  as  the  vertical  shear  in  the 
cross-front  plane,  to  the  forcing  of  the  secondary  circulation.  They  developed  the  theory 
of  describing  temporal  changes  in  the  vertical  and  horizontal  gradients  of  temperature 
and  absolute  momentum  in  the  cross-front  plane  in  terms  of  the  induced  geostrophic  and 
ageostrophic  deformation  field.  It  is  from  this  theory  that  the  Sawyer-Eliassen  diag- 
nostic equation  for  the  secondary  circulation  was  derived.  This  equation  enables  diag- 
nosis of  small-scale  diabatic  and  frictional  processes  as  well  as  those  processes  associated 
with  geostrophic  deformation  forcing  the  secondary  circulation  around  a  front.  Eliassen 
(1962)  demonstrated  that  the  circulation  is  counter-clockwise  (cyclonic)  around  regions 
of  positive  forcing,  and  clockwise  (anti- cyclonic)  around  regions  of  negative  forcing.  The 
circulations  are  concentric  ellipses  in  the  vicinity  of  the  point  source  of  forcing  (Eig.  1). 
oriented  vertically  along  the  front  and  becoming  vertically  or  horizontally  distorted  at 
some  distance  from  the  point  source  as  a  function  of  vorticity  and  inertial  stability 
(Keyser  and  Shapiro  19S6).  The  Sawyer-Eliassen  equation  was  later  expanded  by 
Iloskins  and  Draghici  (1977)  and  Hoskins  et  a!.  ( 1978)  to  account  for  variations  in  the 
along-front  plane,  thus  generalizing  the  theory  of  secondary  circulation  into  three  di- 
mensions, known  as  the  Q-vector  theory.  The  resultant  equation  is  known  as  the 
three-dimensional  (3-D)  form  of  the  Sawyer-Eliassen  equation. 

A.     AGEOSTROPHIC  CIRCULATION  AROUND  A  WARM  FRONT 

Gidel  (197Si  describes  the  processes  of  geostrophic  balance  and  continuity  of  mass 
ilow  which  lead  to  the  secondary  circulation  around  a  warm  front.  The  rising  motion 
within  the  front  coupled  with  the  sinking  motion  on  the  cold  air  side  of  the  front  is  a 
thermodynamically  direct  circulation.  By  continuity  of  mass,  the  How  is  convergent  at 
low  levels  in  the  upward  vertical  motion  region  and  divergent  at  upper  levels.  Transverse 
flow  results  at  lower  and  upper  levels,  completing  the  circulation  cell  in  the  region. 
L'pper-level  westerly  (positive)  momentum  and  low  level  easterly  (negative)  momentum 
are  created,  which  contribute  to  increasing  the  vertical  shear  of  the  along-front  wind 
component.  The  adiabatic  cooling  of  the  rising  warm  air  and  adiabatic  warming  of  the 
sinking  cold  air  contribute  to  decreasing  the  temperature  gradient  across  the  front.    Eigs. 


Fig.  1.  Streamlines  around  n  positive  point  source.  Streamlines  in  the  (a)  m/>-plane 
and  in  the  (b)  jyj-planc.  m  denotes  absolute  momentum.  From  Eliassen 
(1962). 


2  and  3,  from  Gidel  ( 1978).  depict  the  initial  and  advected  potential  temperature  surfaces 
and  isotachs.  and  the  induced  secondary  circulation  around  a  warm  front.  In  the  frontal 
zone  the  wind  field  is  stronger  than  that  required  for  geostrophic  balance,  implying  that 
Coriolis  dominates  the  PGF,  which  sets  up  a  transverse  How  to  the  right  (north).  Im- 
mediately ahead  of  the  front,  the  opposite  is  true,  and  the  resultant  transverse  How  is  to 
the  left  (south).  The  How  is  convergent  directly  beneath  the  front,  and  divergent  behind 
and  ahead  of  the  front.  "I  he  convergence  near  the  leading  edge  of  the  frontal  zone  acts 
as  a  deformation  field  which  enhances  the  potential  temperature  gradient  across  the 
front.  The  circulation  cell  behind  the  front  is  thermodynamicallv  indirect  and  is  the  in- 
direct cell  discussed  by  Eliassen  (1962)  with  respect  to  warm  fronts.  In  the  indirect  cir- 
culation, the  cross-front  thermal  gradient  is  enhanced  by  adiabatic  processes  a'td  the 
along-front  vertical  wind  shear  is  diminished  by  low-level  divergence  and  upper-level 
convereence. 


Idealized  cross  section  through  a  warm  front.  The  initial  position  of  the 
potential  temperature  surfaces  is  indicated  by  the  solid  lines,  and  their  po- 
sition alter  they  have  been  slightly  advectcd  is  indicated  by  dotted  lines. 
The  dashed  lines  are  isotachs  of  the  along-front  wind  field  in  gcostrophic 
balance  with  the  initial  potential  temperature  field.    From  Gidcl  (1978). 
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Fig.  3.  Secondary  circulation  around  a  warm  front.  The  dotted  lines  are  repeated 
from  Fig.  2.  "1  he  dashed  lines  are  the  isotachs  of  the  along-front  wind  field 
after  being  advectcd  from  their  initial  position  shown  in  Fig.  2.  The  solid 
lines  arc  the  isotachs  of  the  along-front  wind  field  u  which  geostrophically 
balances  the  potential  temperature  field,  shown  as  dotted  lines.  The  heavy 
arrows  indicate  horizontal  and  vertical  components  of  the  induced  circu- 
lation.   From  Gidel  (1978). 


The  ageostrophic  flow  around  a  positive  center  of  forcing  is  counter-clockwise, 
which  gives  the  sense  of  the  secondary  direct  circulation  cell  if  it  was  located  ahead  of 
the  warm  front.  Similarly,  a  negative  (clockwise)  center  gives  the  direction  of  flow  of  the 
indirect  cell  if  the  center  was  located  behind  the  front.  Fig.  4a  is  a  depiction  of  the  the- 
oretical location  of  the  positive  and  negative  ellipses  of  forcing  around  a  warm  front. 
A  dipole  effect  is  set  up  between  the  two  centers,  and  the  gradient  of  forcing  from  one 
cell  to  the  other  is  depicted  as  strongest  across  the  vertical  flow  regions,  assuming  the 
downward  and  upward  motions  arc  the  same  magnitude.   Theoretically,  the  stronger  the 


magnitude  and  gradient  of  the  cells,  the  stronger  the  magnitude  offlow  of  the  circulation 
cell.  If  the  updraft  region  is  stronger  than  the  downdraft  region,  one  would  expect  the 
gradients  to  be  as  in  Fig.  4b. 

The  above  discussion  assumes  the  scale  of  geostrophic  and  diabatic  forcing  is  the 
same  as  that  of  the  circulation  cell.  However,  processes  which  contribute  to  enhance- 
ment of  the  circulation  cell  arc  not  on  the  same  scale.  The  Sawyer-Eliassen  equation  is 
a  means  by  which  the  small-scale  processes  of  friction  and  diabatic  heating  can  be  in- 
cluded with  geostrophic  deformation  to  diagnose  the  forcing  of  the  cell.  Total  forcing 
would  not  necessarily  be  a  single  symmetric  center,  but  rather  a  center  distorted  by 
smaller  regions  of  positive  and  negative  forcing  acting  to  enhance  or  retard  the  second- 
ary circulation.  For  example,  one  would  expect  enhanced  regions  of  forcing  at  low  levels 
due  to  geostrophic  and  frictional  confluence,  and  at  mid-levels  due  to  latent  heat  release. 


Fig.  4.  Ageostrophic  forcing  centers  around  a  warm  front.  Ellipses  of  forcing  arc 
for  (a)  an  updraft  and  downdraft  of  equal  strength,  and  (b)  a  stronger  up- 
draft. Note  the  orientation  of  the  front  and  x-  and  y-axes  have  been  re- 
versed from  Figs.  2  and  3. 


B.     SAWYER-ELIASSEN  EQUATION 

Eliassen  (1962)  introduced  the  quantity  absolute  momentum  to  describe  the  second- 
ary circulation  in  terms  of  an  absolute  frame  of  reference,  the  equations  for  which  can 
be  written  as 


mx  =  ug-fy  (1) 

and 

m}  =vg+fx  (2) 

The  thermal  wind  relationships  in  pressure  coordinates  are 


and 


(3) 


M) 


Following  Moskins  and  Draghici  (1977),  and  Keyser  and  Shapiro  (1986),  prognostic 
equations  for  absolute  momentum  and  thermal  advection  can  be  written  as 
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F.  =  — ; — .  and  F.  =  — — .  where  1  is  the  stress  on  a  horizontal  surlacc. 
cp  -        cp 


Q  =  -tj-,  and 


fijr  >«P(-g-). 


By  rearranging  the  terms,  tlie  x  and  y  components  of  the  Sawyer-Eliassen  equation 
can  be  written  as 
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By  taking  the  dot  product  of—  A'  X  yV(?  and  the  vector  equation  made  up  of  the  x 
and  y  components  given  above,  the  magnitude  of  the  forcing  and  circulation  can  be 
obtained.    The  right  hand  side  (forcing)  of  the  equation  becomes 
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The  terms  are  then  normalized  by  dividing  by  |  k  X  ;\0  |,  to  obtain  the  magnitudes 
of  the  components  contributing  to  total  kinematic  forcing  of  the  warm  front.  The 
Savvy  er-Eliassen  equation  is  not  solved;  rather  this  form  enables  a  diagnosis  of  the  con- 
ns to  the  forcing  of  the  circulation  cell  by  each  term.  Centers  of  positive  and 
negative  forcing  are  derived  by  this  diagnosis,  and  represent  centers  of  counter-clockwise 
and  clockwise  circulation  respectively.  The  3-D  form  is  used  in  place  of  the  2-D  version 
because  the  data  set  is  three-dimensional  and  because  the  along-front  gradients  arc  not 
assumed  to  be  in  geostrophic  balance. 

C.     DESCRIPTION  OT  TERMS 

1.     Geostrophic  Confluence 

The  first  term  in  equation  (11)  represents  the  effect  that  confluent  geostrophic 

deformation  lias  in  enhancing  the  horizontal  temperature  gradient  (VpT)  by  compressing 

the  isotherms  across  the  front.  The  strongest  centers  of  confluence  (positive  forcing)  and 

diffluence  should  be  as  in  Fig.  5.    The  low-level  confluence  and  upper-level  dillluence 

represent  a  downward  transport  of  momentum,  which  acts  to  reduce  the  vertical  wind 

shear,  while  at  the  same  time  the  thermal  gradient  is  increasing,    hollowing  Keyser  and 

Shapiro  1 1986),  the  resultant  contributions  to  d(  -~ —  )  /  dt  and  d(y-rr~  !  /  dt  are  »'espec- 

cP   ■  cy 

tiveh  negative  and  positive,  which  is  a  tendency  towards  imbalance  in  the  thermal  wind. 

a  frontogenetic  process.    This  imbalance  is  counteracted  by  the  secondary  direct  circu- 

rced  by  the  confluence,  which  simultaneous!}  increases  the  vertical  wind  shear 


through  the  net  upward  transport  of  positive  momentum,  while  reducing  the  temper- 
ature gradient  by  adiabatic  processes. 
2.     Geostrophic  Shear 

The  second  term  represents  the  effect  that  geostrophic  shear  deformation  has  in 
rotating  the  along-front  temperature  gradient  perpendicular  to  the  front.  If  the  winds 
across  the  front  increase  to  a  maximum  at  the  low  level  easterly  jet  on  the  cold  air  side 
(ahead)  of  a  warm  front,  cyclonic  shear  is  present.  If  the  easterly  wind  aboad  of  the  front 
is  stronger  than  that  behind  the  front,  differential  cyclonic  shear  in  the  presence  of 
along-front  decreases  in  potential  temperature  tends  to  enhance  the  cross-front  temper- 
ature gradient.    Larcer  negative  values  of  u.  ahead  of  the  front  relative  to  weaker  values 

t  ks 
behind  the  front  (— — <  0)  differentially  advect  higher  values  of  potential  temperature 

cy  ' 

at  a  greater  rate  north  of  the  front  relative  to  south  of  the  front.    At  the  same  time,  the 

vertical  wind  shear  (— — )  tends  to  tip  the  ut  isotachs  into  the  vertical  where  the  horizontal 

cp 
shear  is  cyclonic, 'thus  reducing  the  magnitude  ol  the  vertical  wind  shear,  creating  an 

imbalance.   Again,  the  direct  circulation  would  counteract  the  imbalance  created  by  the 

shear  as  with  the  confluence.    One  would  expect  that  the  strongest  forcing  due  to  shear 

would  be  between  the  low-level  easterly  jet  and  the  front. 
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Fig.  5.        Transverse  ageostrophic  circulations  for  an  idealized  front.     Circulations 

o\\  cvg 

characterized  bv  (a)  confluence  (  -r-2-  <  0,  —r—  =  0)  and  (b)  horizontal  shear 

cVj  cvg    '  cy  cp 

(— —  =  0,  — —  >  0).    Dashed  lines,  isotachs  of  ug  (denoted  by  U);  thin  solid 

lines,    isotachs    of  vg   (denoted   by   V);   thick   solid    lines,   streamlines    of 

ageostrophic  circulation.   Adapted  from  Eliassen  (1%2).    From  Keyscr  and 

Shapiro  (-19S6). 


3.     Friction 

The  third  terra  represents  the  convergent  effect  that  friction  has  on  the 
geostrophic  wind.  Winds  at  the  surface  experience  friction,  but  are  maintained  in  part 
by  the  downward  momentum  flux  from  stronger  winds  above  the  surface.  Differential 
downward  momentum  fluxes  act  like  confluence  to  compress  the  temperature  gradient. 
As  the  cyclone  intensifies  the  winds  increase,  which  leads  to  increased  friction  and  fric- 
tional  turning  of  the  winds  toward  the  low.  This  enhances  frictional  convergence.  If  the 
winds  were  to  increase  uniformly  around  the  low,  there  would  be  no  differential  mo- 
mentum flux  across  the  front.  However,  the  wind  gradient  across  a  frontal  zone  is 
greater  than  elswhere  around  the  low.  so  the  maximum  differential  friction  is  in  the  vi- 
cinity of  the  fronts. 

As  discussed  in  Section  I.  favorable  surface  heat  fluxes  also  enhance  frictional 
convergence.  Regions  of  upward  heat  flux  into  the  PBL  from  a  relatively  warm  surface 
can  be  likened  to  warm  air  advection  into  the  region.  This  would  cause  a  veering  of  the 
PBL  wind  wndi  height,  or  an  increased  turning  angle  and  further  convergence  at  the 
surface.  Similarly,  a  downward  heat  flux  from  the  PBL  to  a  cool  surface  would  cause  a 
backing  of  the  thermal  wind  with  height,  or  a  reduced  angle  of  turning  and  reduced 
convergence.  I  or  a  warm  front  in  which  there  arc  downward  heat  fluxes  in  the  cold  air 
sector  ahead  of  the  warm  front  as  well  as  in  the  warm  air  sector  behind  the  front,  one 
can  expect  that  lor  a  given  SST  distribution,  the  air-sea  temperature  difference  would 
be  create!'  in  the  warm  air  region  behind  the  front  than  in  the  cool  air  region  ahead  of 
the  front.  '1  he  downward  heat  flux  would  be  greater  behind  the  front  than  ahead,  with 
tiie  backing  of  the  thermal  wind  stronger  behind  the  front  than  ahead.  This  creates  a 
cross  the  front,  with  stronger  convergence  ahead  of  the  front.  "1  litis,  the 
frictional  convergence  is  composed  of  the  geostrophic  component  and  a  heat  flux  com- 
ponent, fig.  o  depicts  the  geostrophic  frictional  forcing  as  a  result  of  the  differential 
downward  transfer  of  momentum  from  the  winds  above  the  PBL.  and  the  heat  flux 
component  of  frictional  forcing  within  the  PBL  due  to  surface  heat  fluxes.  The  depiction 
is  for  the  case  of  stronger  downward  momentum  fluxes  behind  the  front  than  ahead,  and 
for  stronger  downward  heat  fluxes  at  the  surface  behind  the  front.  This  would  result  in 
a  reinforcement  of  frictional  confluence.  The  greatest  differential  would  be  across  the 
warm  front,  and  the  strongest  frictional  confluence  would  be  located  where  the  strongest 
differential  occurs.  Realistically  it  is  possible  to  have  stronger  downward  momentum 
fluxes  ahead  of  the  warm  front,  which  would  oppose  the  heat  flux  component  depicted 
above,  resulting  in  neelicible  frictional  forcing. 
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Fig.  6.  Centers  of  frictional  forcing.  Heavy  lines  above  PI? L  denote  secondary 
circulation  cells,  with  vertical  arrows  in  the  PBL  denoting  downward  mo- 
mentum fluxes  (solid)  and  downward  heat  fluxes  (curved).  Heavy  arrows 
indicate  stronger  downward  momentum  fluxes.  Dashed  lines  depict  warm 
front. 


4.     Diabatic  Heating 

The  fourth  term  represents  the  effect  diabatic  heating  has  on  the  horizontal 
temperature  gradient.  In  the  boundary  layer,  differential  heat  fluxes  act  to  relax  the 
temperature  gradient,  and  arc  strongly  dependent  upon  the  surface  temperature  distrib- 
ution, as  discussed  in  Section  I.  Surface  heat  and  moisture  fluxes  in  the  PBL  can  be 
turbulently  mixed  out  of  the  PBL,  providing  a  source  of  heat  and  moisture  to  the  mid- 
troposphere  in  addition  to  the  heat  and  moisture  advectcd  into  the  region.  At  mid- 
levels,  latent  heat  release  within  and  behind  the  warm  front,  and  evaporative  cooling 
ahead  of  the  front  enhance  the  temperature  gradient.  One  should  expect  centers  of 
positive  and  negative  diabatic  forcing  to  be  located  as  in  fig.  7.  I  he  low  level  forcing 
opposes  the  secondary  circulation  at  low  levels  and  the  mid  level  forcing  enhances  the 
circulation. 


Fig.  7.        Centers  of  diabalic  forcing.    Forcing  in  the  PHI.  opposes  the  secondary 

circulation.    In  the  mid-troposphere  it  enhances  the  circulation. 


III.     DESCRIPTION  OF  DATA  SET 

A.     NUMERICAL  MODEL 

The  numerical  simulations  for  the  three  case  studies  used  in  this  study  were  run  for 
4S  hours  using  the  Pennsylvania  State  University  National  Center  for  Atmospheric  Re- 
search (PSU  NCAR)  mesoscale  model  version  4  (MM4),  which  is  described  in  detail  by 
Anthes  et  al.  (19S7).  The  usual  meteorological  data  plus  the  boundary  layer  fluxes  were 
output  every  3  hours.  The  version  of  the  model  used  for  this  study  utilized  14  vertical 
layers  from  the  surface  to  100  mb,  periodic  east-west  boundaries,  a  ft  -  plane  variation 
of  the  Coriolis  parameter,  and  the  high  resolution  boundary  layer  model.  The  boundary 
layer  occupies  five  of  the  14  layers,  with  the  lowest  layer  within  the  surface  layer. 

Differential  friction  can  result  from  the  boundary  layer  structure  and  processes  in 
the  model  since  the  model  vertical  heat,  moisture,  and  momentum  fluxes  use  a  vertical 
diffusion  coefficient  that  depends  on  the  local  bulk  Richardson  number  (Rb).  The 
boundary  conditions  are  considered  stable  when  the  bulk  Richardson  number  (Rb)  is 
greater  than  zero  and  ZHjL  is  large.  ZH  i^  the  height  of  the  surface  layer  and  L  is  the 
Monin-Obukhov  length.  For  very  stable  conditions  (#4>0.2  and  ZH\L  unrestricted), 
no  turbulent  fluxes  occur,  for  stable  conditions  (0.0  <  RD  <  0.2  and  ZHjL  unrestricted). 
damped  mechanical  turbulent  fluxes  occur,  and  for  marginally  unstable  or  neutral  con- 
ditions (Rb<  0.0  and  ZHjL  <  1.5),  both  mechanical  and  buoyant  turbulent  fluxes  occur. 
Instable  boundary  conditions  apply  when  both  Rb  <  0.0  and  ZH\L>  1.5,  and  under 
these  conditions  the  vertical  fluxes  are  computed  using  a  convective  plume  model.  In- 
cluded in  the  fluxes  is  a  downward  entrainment  flux  at  the  top  of  the  mixed  layer  equal 
to  20%  ol'  the  surface  flux. 

Cloud  formation  and  precipitation  are  parameterized  through  both  non-convective 
and  convective  processes.  Stable  lifting  that  results  in  relative  humidities  greater  than 
100%  is  treated  as  non-convective  precipitation,  and  the  latent  heat  is  added  directh  to 
the  thermodynamic  equation.  Unstable  ascent  or  convection  is  parameterized  using  a 
modified  kuo  (1974)  scheme  (Anthes  1977)  which  redistributes  the  moisture  and  latent 
heat  vertically  through  prescribed  profiles.  The  criteria  for  convection  are  total  column 
moisture  convergence  greater  than  0.01  g  m  :  ^  '  and  the  predicted  cloud  height  greater 
than  700  mb. 


Diabatic  heating  used  in  the  Sawyer-Eliassen  diagnostic  equation  is  computed  as  the 
residual  of  the  thermodynamic  equation 

<M)  cT  r'V 

(J  =  -V  =  C (  -7^-  +  J" .  Vn T  +  v)  -4-  )  +  ZG) 


0=r(-^)exp-^- 

The  time  tendency  was  computed  as  a  six-hour  time  difference  and  introduces  some  error 
into  the  diabatic  heating  values.  However,  the  diabatic  heating  data  does  correspond 
with  the  rainfall  regions,  giving  at  least  qualitative  confidence  in  their  distribution.  Be- 
cause frictional  forcing  is  computed  as  a  second  order  vertical  differential,  forcing  at  the 
surface  is  represented  by  forcing  at  a  level  above  the  surface  in  the  middle  of  the  PBL. 

B.     INITIAL  CONDITIONS 

The  initial  sea-level  pressure  (SLP),  surface  equivalent  potential  temperature,  and 

observed  wind  field  for  the  three  cases  are  shown  in  Fig.  8.    The  initial  baroclinic  insta- 

nditions  are  the  same  in  each  case,  with  the  only  differences  due  to  the  respective 

:c  heat  and  moisture  flux  distributions.  The  initial  model  cyclone  differs  from  an 
actual  atmospheric  cyclone  most  significantly  in  its  lack  of  an  initial  vertical  circulation, 
its  moisture  structure  being  too  dry  compared  with  typical  cloud  patterns  of  actual 
cyclones,  and  a  slight  upshear  tilt  in  the  thermal  wave.  All  but  the  moisture  distribution 
difference  are  eliminated  in  the  first  six  hours  of  each  model  simulation.  The  moisture 
distribution  takes  IS  hours  to  resemble  the  typical  cloud  pattern  of  an  actual  cyclone, 
which  also  gives  less  certainty  to  the  diabatic  heating  contributions  to  the  model 
cyclogenesis  and  frontogenesis. 

1  he  control  experiment  (Case  1)  uses  a  zonal  SST  distribution  (Fig.  9)  and  includes 
surface  heat  and  moisture  fluxes  throughout  the  simulation.  With  this  SST  distribution 
the  mean  baroclinicity  in  the  atmosphere  is  reinforced,  and  the  pattern  of  surface  heating 
introduces  temperature  tendencies  that  counteract  the  initial  thermal  advection  (Nuss 
and  Anthes.  1987).  The  initial  surface  flux  distribution  is  characterized  by  upward  heat 
and  moisture  fluxes  west  of  the  surface  low  and  downward  fluxes  east  of  the  surface  low. 
The  second  experiment  (Case  2)  uses  the  same  zonal  SST  distribution  but  insulates  the 
atmosphere  from  the  sea  surface  by  specifying  a  zero  heat  and  moisture  flux  condition 


Fig.  8.  Initial  conditions.  SI  P  in  nib  (solid).  0t  in  K  (dashed),  and  observed  wind 
vectors  in  in  s  ]  lor  each  case.  Length  of  wind  vector  denotes  intensity, 
with  the  maximum  length  representing  50  in  s-1. 


throughout  the  simulation.  Because  the  only  source  of  moisture  in  the  model  is  from  the 
surface  moisture  flux,  this  case  has  no  moisture  directly  transported  to  mid  and  upper 
levels  as  in  the  other  two  cases,  and  relies  on  the  cloud  formation  and  precipitation 
parameterization  applied  to  rising  warm  air  cooled  through  adiabatic  processes.  The 
downward  momentum  flux  is  not  specifically  modified  from  Case  1,  except  indirectly 
through  stratification  effects  that  result  from  no  surface  heating.  With  no  surface  lluxcs, 
the  total  frictional  convergence  is  reduced  to  that  caused  by  differential  downward  mo- 
mentum fluxes,  and  there  is  no  thermal  damping  effect  in  the  boundary  layer. 
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Fig.  9.  Case  1  and  Case  2  SSI  distribution.  SSI  in  °C (solid),  with  the  initial  total 
(latent  plus  sensible)  licit  flux  distribution  in  \V  m  -\  and  the  cyclone  track 
and  position  at  six-hour  intervals  from  00  h  to  48  h  indicated  for  Case  1. 


The  final  experiment,  Case  3,  uses  a  sinusoidal  SST  distribution  and  includes  surface 
heat  and  moisture  fluxes  as  in  <  asc  1  (Fig.  10).  The  sinusoidal  SST  distribution  in  this 
case  is  in  phase  with  the  thermal  structure  of  the  cyclone,  thus  reinforcing  the  low-level 
baroclinicity  instead  of  the  large  scale  mean  baroclinicity.  This  distribution  results  in 
large  upward  total  heat  fluxes  to  the  northeast  of  the  surface  low  as  well  as  to  the  west. 
In  this  case  the  magnitude  of  the  heat  flux  is  largest  northeast  of  the  surface  low,  en- 
hancing the  frictional  convergence  northeast  of  the  surface  low  compared  with  Case  1. 


Fig.   10.        As  in  Fig.  9,  except  for  Case  3. 


C.     DEVELOPMENT  On  I  IF  CYCLONE 

The  differences  in  structure  and  dynamics  of  each  cyclone  throughout  the  simu- 
lations represent  the  direct  and  indirect  effects  of  the  surface  ilux  distribution  unique  to 
each  case,  since  the  initial  baroclinic  instability  factors  arc  identical.  The  intensity  of 
cyclogenesis  in  each  case  can  be  described  in  part  by  the  sea-level  pressuie  (SEP) 
deepening  rate.  Fig.  11  is  a  plot  of  the  lowest  closed  contour  of  SLP  around  the  center 
versus  time  for  each  case.  By  the  Sanders  and  Gyakum  (1(>S())  definition  of  explosive 
development  (1  nib  /r'  for  24  h  or  more),  all  three  cyclones  develop  explosively,  with  the 
a\erage  deepening  rate  1.40  nib  //  '  for  Case  1,  1.89  nib  h  ]  for  Case  2,  and  1.65  nib  h  ' 
for  Case  3  (Nuss  and  Anthes  1987). 

In  Case  1,  the  cyclone  tracks  northeast  towards  colder  water,  initially  over  n  region 
of  increasing  downward  heat  fluxes,  and  then  towards  the  region  of  increasing  upward 
heat  fluxes  in  the  top  right  quadrant  of  Fig.  9.  As  the  cyclone  moves  towards  this  re- 
gion, the  differential  in  static  stabilitv  across  the  warm  front  is  enhanced.     Increased 


Fig.  11.  SLP  (ml))  Aeisus  time  for  each  case.  Note  that  the  deepening  rates  are 
the  same  lor  each  case  until  12  h,  and  the  deepening  rates  of  Cases  2  and 
3  arc  equivalent  until  30  h. 


upward  heat  (luxes  ahead  of  the  front  reduce  the  static  stability  in  this  region  compared 
with  the  comparatively  small  increase  in  static  stability  behind  the  front  in  the  region 
of  downward  heat  fluxes.  Dilfciential  static  stability  in  the  PBL  increases  fiictionai 
convergence,  however  it  is  small  relative  to  the  thermal  damping  effects  of  the  heat  flux 
distribution,  which  acts  to  slow  fiontogcncnsis  by  warming  the  cold  air  sector  ahead  of 
the  front  and  cooling  the  warm  air  sector  behind  the  front.  'J  his  prevents  the  formation 
of  concentrated  temperature  gradients  at  the  surface  and  associated  strong  temperature 
advection  (Nuss  and  Anthcs,  1VS7). 

The  removal  of  surface  heat  and  moisture  fluxes  in  Case  2  results  in  an  increased 
rate  of  cyclogcnesis  compared  with  Case  1.  In  the  absence  of  surface  fluxes,  differential 
static  stabilty  across  the  front  is  reduced,  so  frictional  convergence  due  to  differential 
downward  momentum  fluxes  is  no  longer  enhanced.  However,  the  frontolytical  thermal 
effect  of  the  Case  1  surface  heat  fluxes  is  removed,  and  indeed.  Nuss  and  Anthes  (19S7) 
showed  that  the  low-level  thermal  cold  and  warm  front  advection  increased  in  this  case 


compared  with  Case  1.  and  the  temperature  gradients  were  stronger,  indicating  stronger 
cold  and  warm  fronts. 

In  Case  3  the  cyclone  tracks  northeast  across  the  SST  gradient,  moving  slightly  to- 
wards warmer  water  until  the  24  hour  point  (24  h)  when  the  cyclone  starts  tracking  to- 
wards colder  water.  This  track  results  in  an  initial  reduction  in  downward  fluxes  east 
of  the  surface  low  until  24  hours,  when  the  downward  fluxes  begin  to  increase  as  the 
cyclone  changes  direction.  The  surface  flux  distribution  of  the  first  24  hours  favors  a 
maintenance  of  the  baroclinicity  of  the  cyclone.  The  distribution  of  strong  upward 
fluxes  north  of  the  warm  front  results  in  reduced  low  level  static  stability  and  stronger 
heating  of  the  PBL  in  this  region,  which  enhances  frictional  convergence  northeast  of  the 
surface  low.  The  frictional  convergence  in  this  case  is  stronger  than  in  Case  1.  where  the 
upward  heat  fluxes  northeast  of  the  surface  low  were  not  as  strong.  Another  significant 
difference  between  the  two  cases  is  that  the  moisture  flux  to  the  northeast  of  the  warm 
front  in  Case  3  is  greater  than  200  W  m~2  compared  with  near  zero  in  the  control  case 
(Nuss,  1989).  and  there  is  an  upward  moisture  flux  throughout  most  of  the  warm  sector 
in  Case  3  compared  with  near  zero  fluxes  in  Case  I. 

As  discussed  by  Nuss  (1989),  convergence  near  the  cyclone  center  and  warm  front 
p!  iy  a  significant  role  in  the  vertical  circulation  of  the  cyclone.  Because  of  the  heat  flux 
distribution,  stronger  wind  stress  results  to  the  northeast  of  the  surface  low  in  Case  3 
compared  with  Cases  1  and  2.  The  curl  of  the  surface  stress  is  proportional  to  the  fric- 
tional PBL  convergence  (Fleagle  and  Nuss.  19S5).  With  increased  convergence  into  the 
surface  low  and  warm  front  in  Case  3  as  a  result  of  the  stronger  heating,  the  PBL 
ageostrophic  flow  toward  the  front  and  low  is  strongest  in  this  case.  This  results  in  a 
stronger  low-level  transverse  component  of  the  secondary  circulation,  which  in  turn  en- 
hances the  intensity  of  the  cold  conveyor  belt  flowing  towards  the  cyclone  ahead  of  the 
warm  front,  which  transports  boundary  layer  air  that  has  been  heated  and  moistened 
by  the  surface  heat  fluxes  toward  the  cyclone  center.  The  result  is  increased  vertical 
transport  of  heat  and  moisture  out  of  the  boundary  layer,  contributing  to  stronger 
midtropospheric  heating  and  greater  intensification  of  the  cyclone  than  in  Case  1.  In 
both  Cases  1  and  3,  the  strongest  ascent  was  found  to  be  along  a  line  extending  from  the 
cyclone  center  southeast  and  along  the  warm  front.  In  Case  3.  both  the  frictional  con- 
vergence and  latent  heat  release  are  offset  by  the  frontolytical  effect  of  the  surface  heat 
flux  distribution  after  24  h.  It  appears  that  the  reduction  in  the  deepening  rate  of  the 
cyclone  can  be  related  to  the  removal  of  the  favorable  conditions  of  the  first  24  hours. 
It  remains  to  be  shown  whether  the  stronger  frictional  convergence  in  Case  3  is  due  to 


PBL  processes  or  stronger  geostrophic  forcing  that  results  from  other  indirect  effects  of 
the  surface  fluxes. 

The  role  of  surface  heat  and  moisture  fluxes  in  oceanic  cyclogenesis  is  complex.  The 
frontogenctical  effects  of  surface  heat  fluxes  on  the  PBL  stratification  and  frictional 
convergence  are  opposed  by  the  direct  frontolytical  effects  on  the  thermal  gradient.  On 
the  other  hand,  turbulent  mixing  can  transport  heat  and  moisture  out  of  the  boundary 
layer  into  the  troposphere  where  it  can  enhance  the  temperature  gradient  through  the 
release  of  latent  heat.  This  feeds  back  to  the  cyclone,  since  the  enhanced  thermal  gra- 
dient enhances  the  secondary  circulation  and  intensity  of  the  cyclone,  which  in  turn  en- 
hances the  winds  and  the  geostrophic  frictional  convergence.  It  can  be  seen  that  the 
advance  of  the  cyclone  into  a  new  distribution  of  differential  heat  and  moisture  fluxes 
has  a  critical  impact  on  the  intensity  of  cyclogenesis.  Case  1  allows  an  analysis  of 
boundary  layer  processes  in  the  presence  of  a  continuously  changing  cross-front  SST 
distribution.  Case  2  allows  analysis  of  frontogenesis  and  cyclogenesis  in  the  absence  of 
m  isture  fluxes,  and  Case  ?  allows  analysis  of  these  processes  in  the 
presence  of  an  initially  uniform  SST  gradient  followed  by  a  changing  gradient  alter  24 

D.  TOTAL  FORCING  AT  THL  SURFACE 

.  itensity  of  the  surface  warm  front  can  be  represented  by  the  total  geostrophic 
and  agcostrophic  forcing  at  the  surface.  I  igs.  12  and  13  show  the  total  forcing  at  the 
surface  at  24  h  lor  Cases  1  and  2.  and  Fig.  14  shows  the  same  for  Case  3.  Recalling  the 
lireci  '   around  vertically-oriented  centers  of  positive  and  negative  forcing  as 

discussed  in  Section  II.  one  can  determine  the  sense  of  flow7  in  the  vertical  plane  from 
the  sign  of  the  forcing  at  the  surface.  In  each  case  there  is  a  region  of  slightly  negative 
forcing  ahead  of  the  warm  front,  which  implies  clockwise  circulation  in  the  vertical  plane 
and  consequently  flow  towards  the  northeast  in  the  horizontal  plane  of  the  surface. 
Similarly,  the  regions  of  positive  forcing  imply  counterclockwise  circulation  in  the  verti- 
cal plane,  and  consequently  surface  flow  towards  the  southwest.  The  vertical  plane  of 
flow  is  assumed  to  be  perpendicular  to  the  axis  of  horizontal  forcing,  and  the  sign  of  the 
forcing  at  the  surface  is  assumed  to  extend  vertically  above  that  region. 

The  boundaries  between  the  positive  and  negative  regions  of  forcing  are  of  particular 

importance  since  they  represent  the  transition  from  counterclockwise  to  clockwise  How. 

and  hence  are  regions  of  upward  or  downward  vertical  motion.    For  example,  in  Case 

undan  between  the  negative  region  to  the  northeast  of  the  front  and  the  posi- 


tive  region  within  the  surface  front  is  a  downward  vertical  motion  region.  Similarly,  the 
boundary  between  this  positive  region  and  the  negative  region  south  of  the  surface  low 
is  an  upward  vertical  motion  region.  This  combination  of  downward  motion  ahead  of 
the  front  and  upward  motion  behind  the  front  supports  the  direct  circulation  and  is 
frontogenctical  in  that  the  forcing  enhances  the  temperature  gradient.  The  strongest 
surface  flow  is  to  the  southwest,  perpendicular  to  the  axis  of  forcing,  which  is  close  to 
the  position  of  the  warm  front. 

The  total  forcing  is  strongest  in  Case  2  and  weakest  in  Case  1.  which  supports  the 
results  of  Nuss  ( 1989).  who  found  that  the  Case  2  warm  front  was  the  strongest.  Yet  this 
case  did  not  have  the  strongest  vertical  circulation  or  cold  conveyor  belt,  nor  did  it  have 
the  most  rainfall,  although  it  had  the  most  intense  cyclone.  The  structure  of  the  vertical 
circulation  and  forcing  above  the  surface  is  not  necessarily  the  same  as  that  of  the  sur- 
face. A  clearer  picture  can  be  obtained  from  an  analysis  oC  the  cross-sections  of  total 
forcing  for  each  case.  Of  more  importance  is  how  each  component  of  forcing  is  con- 
tributing to  the  total  forcing  at  all  levels,  and  what  effect  the  surface  fluxes  have  on  these 
components. 


Fig.  12.  Total  forcing  at  the  surface  at  2-4  h,  Case  1.  Geopotential  height  in  m 
(solid),  total  forcing  in  K  100  knrx  day1  (dashed),  and  observed  wind 
vectors  in  m  s  '.  Wind  vectors  vary  in  length  according  to  speed,  with  the 
maximum  length  50  in  s  '.  Contour  interval  of  forcing  is  5.  Heavy  lines 
denote  zero  contours. 
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Tig.   13.       As  in  Fig.  12.  except  for  Case  2.   Contour  interval  is  10. 


Tig.   14.       As  in  Fig.  12.  except  for  Case  3.   Contour  interval  is  1C 


IV.     RESULTS 

A.  OVERVIEW 

A  scries  of  surface  plots  of  total  frontogcnctical  forcing  through  the  cyclone  domain 
were  taken  at  six-hour  intervals  to  analyze  the  evolution  of  total  forcing.     Figs.   12 

through  14  are  representative  of  the  general  pattern  and  relative  intensities  of  forcing 
of  the  cold  and  warm  fronts  observed  throughout  the  simulations.  Although  there  is 
positive  forcing  of  the  warm  front  in  each  case,  total  forcing  at  the  surface  in  Cases  1 
and  3  is  strongest  in  the  cold  air  sector  behind  the  cold  front,  whereas  in  Case  2  total 
forcing  is  strongest  in  the  vicinity  of  the  warm  front.  This  suggests  that  processes  unique 
to  Cases  1  and  3  are  dominating  total  forcing  at  the  surface. 

B.  TOTAL  FORCING  OF  THE  WARM  FRONT 

1.     Surface 

As  seen  in  Figs.  12  through  14,  the  pattern  of  forcing  of  the  warm  front  in  Cases 
2  and  3  is  elliptical,  whereas  a  distinguishable  ellipse  of  forcing  does  not  become  evident 
until  30  h  in  Case  1.  Total  forcing  reaches  a  maximum  at  30  h  in  Case  3.  and  36  h  in 
Cases  1  and  2.  Fig.  15  is  a  time  series  of  the  total  forcing  of  the  warm  front,  measured 
as  the  maximum  contour  value  at  the  approximate  mid  point  of  the  surface  warm  front. 
This  is  in  the  location  of  maximum  forcing,  and  also  avoids  the  occluded  front.  Total 
forcing  reaches  a  maximum  of  10  K  1(H)  km  l  day-1  in  Case  1,  40  K  100  km~]  day1  in 
Case  3.  and  100  K  100  km'1  day-1  in  Case  2.  Comparing  the  SLP  deepening  rates  from 
Fig.  1 1  with  the  rates  of  growth  and  relative  intensities  of  total  forcing  from  Fig.  15  for 
each  case,  one  can  see  that  the  strongest  frontogcnctical  forcing  and  deepest  cyclone  are 
both  in  Case  2,  and  the  weakest  total  forcing  and  slowest  deepening  rate  are  in  Case  1. 
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Evolution  of  total  forcing  of  the  surface  warm  front. 


2.     Cross  Sections 

Cross  sections  of  total  forcing  of  the  warm  front  were  taken  at  six-hour  intervals 
for  each  case,  oriented  perpendicular  to  the  warm  front  and  bisecting  the  region  oC 
maximum  total  forcing  of  the  warm  front  at  the  surface  for  that  time.  The  cross  sections 
depict  more  clearly  the  location  of  total  forcing  in  the  vertical  plane  with  respect  to  the 
vertical  circulation  field  and  surface  warm  front.  By  12  hours  all  three  cases  have  broad 
regions  of  positive  total  forcing  in  the  cold  air  sector  ahead  of  the  warm  front,  and  gen- 
erally negative  forcing  in  the  warm  air  sector  behind  the  warm  front.  The  counter- 
clockwise circulation  induced  by  the  positive  forcing  ahead  of  the  front,  and  the 
clockwise  circulation  induced  by  the  negative  forcing  behind  the  front,  set  up  the  direct 
circulation  cell  and  upward  vertical  motion  region  in  the  warm  front,  as  suggested  in 
Section  II  and  as  indicated  in  Figs.  16  through  18  by  the  vectors  of  agcostrophic  vertical 
motion.  Figs.  16  and  17  show  the  cross  sections  of  total  forcing  at  24  hours  for  Cases 
1  and  2.  and  Fig.  IS  shows  the  same  for  Case  3.  The  orientation  from  left  to  right  in  the 
cross-sections  is  northeast  to  southwest,  with  the  cold  air  sector  ahead  of  the  front  to 
the  left  and  the  warm  air  sector  behind  the  front  to  the  right. 

The  initial  model  cyclone  does  not  have  an  initial  vertical  circulation  field;  rather 
it  evolves  over  the  course  of  the  simulation  as  frontogenctical  processes  occur  in  the 
model.      The   presence   of  a   distinct   updraft   region,   as   depicted   by   the   vectors   of 


ageostrophic  flow  around  the  warm  front,  is  evident  by  the  12  hour  point  in  Case  3,  with 
a  less  distinct  updraft  region  in  Case  1  at  this  time.  Case  2  docs  not  develop  a  distin- 
guishable updraft  region  until  IS  hours.  One  can  get  a  sense  of  the  intensity  of  the  cold 
conveyor  belt  (CCB)  from  the  magnitude  of  the  ageostrophic  vectors  circulating  around 
the  direct  cell  in  each  case.   The  strongest  CCB  is  in  Case  3,  while  the  weakest  is  in  Case 
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Fig.  16.  24-h  cross  section  of  total  ageostrophic  forcing,  Case  1.  forcing  is  in  K 
100  knr1  day1  (solid).  0,  in  K  (dashed),  and  ageostrophic  circulation  vec- 
tors. Ageostrophic  circulation  vectors  vary  in  length  according  to  speed, 
with  the  maximum  length  corresponding  to  50  m  s  '  in  the  horizontal  and 
50  /.ib  s-1  in  the  vertical.  Heavy  lines  denote  zero  contours.  Contour  in- 
terval of  forcing  is  5. 
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Fig.   17.        As  in  Fig.  16,  except  for  Case  2.    Contour  interval  is  5. 
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fig.   18.       As  in  Fig.  16,  except  for  Case  3.   Contour  interval  is  10. 


By  the  18  hour  point  Cases  1  and  3  have  distinguishable  ellipses  of  positive  total 
forcing  extending  vertically  within  and  ahead  of  the  updraft  region,  whereas  in  Case  2  a 
similar  region  is  not  apparent  until  24  hours.  Cores  of  positive  forcing  in  the  boundary 
layer  are  first  distinguishable  in  Case  1  at  30  h,  and  at  24  h  in  Cases  2  and  3,  after  which 
the  boundary  layer  remains  the  level  at  which  total  forcing  is  the  strongest.  Of  all  the 
cases,  Case  2  develops  the  strongest  boundary  layer  total  forcing,  reaching  a  maximum 
of  150  K  100  knr1  dqy~]  at  36  h.  The  relative  maxima  and  times  of  occurrence  in  Cases 
1  and  3  are  35  K  100  km~l  day  '  at  36  h  and  60  K  100  knrx  day  ]  at  30  h  respectively. 
Tig.  19  shows  the  evolution  of  total  forcing  in  the  BBL  for  each  case,  measured  from  the 
respective  cross  sections  as  the  maximum  positive  value  in  the  core  of  positive  forcing 
near  the  updraft  region,  typically  at  about  950  mb.  Although  this  is  a  crude  measure 
of  the  frontogenesis,  it  gives  at  least  qualitative  evidence  of  how  the  forcing  evolves  in 
each  case.  Case  3  has  the  strongest  initial  forcing,  but  is  the  case  in  which  the  PBL 
forcing  first  declines.    Between  IS  h  and  24  h  the  total  forcing  in  Cases  2  and  3  is  the 


same,  after  which  the  rate  of  growth  begins  to  decline  in  Case  3  while  the  growth  rate 
in  Case  2  continues  to  increase. 
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Fi».   19.        Evolution  of  total  forcing  in  the  PBL. 


At  mid  levels.  Case  1  develops  the  strongest  total  forcing,  reaching  35  K  100 
knrx  dayx  at  42  h.  Case  2  develops  a  maximum  of  30  K  100  km  '  day-\  and  Case  3 
develops  a  maximum  of  2c  K  1«>*»  km  '  </<n~'  at  24  h.  Fig.  2<>  is  a  summary  of  the  total 
forcing  at  mid  levels,  measured  in  the  direct  circulation  cell  as  the  maximum  contour  of 
positive  forcing  at  600  mb. 
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Fig.  20.        Evolution  of  total  forcing  at  mid  levels. 


C.     TOTAL  FORCING  AT  2-4  I  lOURS 

For  continuity  with  Nuss  (19S9),  the  24-hour  time  was  chosen  for  detailed  analysis. 
This  time  is  representative  of  the  most  rapid  phase  of  development,  and  is  just  prior  to 
formation  of  the  occluded  front.  By  this  time  in  the  evolution,  the  effect  of  the  fluxes 
will  have  influenced  all  the  terms  to  some  degree.  In  all  three  cases  there  is  a  large  region 
of  positive  forcing  ahead  of  the  warm  front,  and  generally  negative  forcing  behind  the 
front.  The  gradient  of  total  forcing  at  mid  levels  is  strongest  in  the  updraft  region  where 
the  vertical  vectors  are  the  strongest.  The  positive  forcing  in  the  direct  circulation  cell 
dominates  the  negative  forcing,  suggesting  that  the  direct  cell  is  the  dominant  cell  in  the 
secondary  circulation  of  the  warm  front. 

The  patterns  of  forcing  of  the  direct  circulation  cell  in  each  case,  as  seen  in  Figs. 
16  through  18,  arc  not  symmetrical  ellipses  as  discussed  in  Section  II,  particularly  in 
Cases  2  and  3  where  there  is  a  cell  of  strong  positive  forcing  in  the  boundary  layer  be- 
neath the  updraft  region.  Case  1  develops  a  similar  region  by  30  hours.  There  arc  dis- 
tinct dilferenccs  in  the  forcing  at  each  level  in  the  three  cases,  and  as  discussed  in  the 
previous  section,  must  be  a  result  of  the  diflcrcnccs  in  the  surface  (luxes  in  ea;  h  case. 
As  described  by  Nuss  (1989),  different  surface  heat  and  moisture  flux  distributions  re- 
sulted in  differing  rates  of  cyclogcncsis,  frontogenesis.  and  intensity  of  the  cold  conveyor 
belt.     There  is  certainly  evidence  of  a  relationship  between  total  forcing  of  the  warm 


front,  especially  at  the  surface,  and  the  cyclone  deepening  rate.  The  important  question 
is  how  these  flux  distributions,  which  differed  primarily  near  the  warm  front,  modified 
the  warm  frontal  temperature  gradient  and  vertical  motion  in  each  case,  as  indicated  by 
differences  in  total  forcing. 

D.     CONTRIBUTIONS  BY  EACH  COMPONENT  OF  FORCING 

1.  Overview 

The  total  forcing  in  each  case  is  strongest  in  the  boundary  layer,  with  confluence 
and  shear  the  dominant  contributors  to  total  forcing  at  this  level.  Confluent  forcing  is 
generally  concentrated  beneath  the  updraft  region,  with  shear  extending  from  this  lo- 
cation into  the  direct  circulation  cell  at  low  levels.  Diabatic  heating  is  generally  negative 
at  this  level  in  all  three  cases,  concentrated  beneath  the  updraft  region  but  extending 
ahead  of  the  warm  front  as  with  shear.  Frictional  forcing  is  a  small  component  of  total 
forcing  of  the  warm  front.  The  forcing  due  to  diabatic  heating  is  the  dominant  compo- 
nent of  total  forcing  at  mid  and  upper  levels,  and  generally  located  within  and  ahead  of 
the  updraft  region.  Centers  of  mid-level  confluence  and  shear  are  located  within  the  di- 
rect cell  ahead  of  the  warm  front.  Interesting  departures  from  these  general  results  were 
found  in  each  case,  as  described  below. 

2.  Case  1 

The  24-h  sensible  and  latent  heat  flux  distributions  for  Case  1  are  shown  in  Fig. 
21.  with  the  positions  of  the  front  and  cross-section  added  for  reference.  By  this  time 
the  initial  1009  mb  low  has  deepened  to  a  986  nib  low.  with  strong  positive  total  (latent 
phis  sensible)  heat  fluxes  west  of  the  low  and  downward  fluxes  east  of  the  low.  This 
distribution  prevents  the  formation  of  concentrated  surface  temperature  gradients  and 
strong  temperature  advection,  since  the  upward  fluxes  aie  warming  the  cold  air  region 
and  the  downward  fluxes  arc  cooling  the  warm  air  region.  However  this  effect  extends 
only  to  the  top  of  the  boundary  layer  (Nuss  and  Anthes,  1987). 


Fig.  21.  2-4-h  heat  flux  distributions  for  Case  1.  (a)  sensible  and  (b)  latent  heat 
flux  distributions  in  W  m  2.  Dashed  lines  indicate  downward  or  negative 
fluxes.  Fronts  added  for  reference.    From  Nuss,  1989. 


a.     Surface 

The  frictional  convergence  and  diabatic  heating  components  of  forcing  for 
the  surface  are  shown  in  I  igs.  22  and  23.  Frictional  forcing  is  negative  in  the  region  of 
upward  heat  fluxes  northeast  of  the  surface  low,  but  also  slightly  negative  in  the  region 
of  downward  heat  (luxes  behind  the  warm  front.  However,  the  frictional  forcing  is 
strongly  positive  in  the  cold  air  sector  behind  the  cold  front  where  the  upward  heat  flux 
is  the  strongest.  1  he  inconsistency  of  the  frictional  forcing  relative  to  the  sign  ol  the 
heat  fluxes  suggests  that  the  momentum  fluxes  in  the  model  depend  on  other  factors  that 
are  not  determined  simply  from  the  fluxes.  This  may  be  particularly  important  when  the 
boundary  layer  is  close  to  neutral  stratification. 

The  sense  of  the  ageostrophic  flow  induced  in  the  region  of  negative  forcing 
northeast  of  the  warm  front  is  for  surface  flow  to  the  northeast,  with  corresponding 
downward  motion  at  the  boundary  between  this  region  and  the  region  of  positive  forcing 
closer  to  the  front.    Similarly,  the  induced  flow  in  this  region  of  positive  forcing  is  to  the 


southwest,  with  downward  flow  along  the  boundary  with  the  negative  region  to  the 
northeast  closer  to  the  warm  front.  The  boundary  between  this  positive  forcing  region 
and  the  region  of  negative  frictional  forcing  associated  with  the  low  and  warm  front  is 
upward.  The  combined  flow  would  tend  to  support  the  direct  circulation  cell  as  dis- 
cussed in  the  previous  section.  The  overall  pattern  of  frictional  forcing  resembles  the 
pattern  of  total  heat  fluxes.  Over  the  course  of  the  simulation  the  band  of  positive  fric- 
tional forcing  expands  in  size,  although  the  magnitude  does  not  substantially  change. 
This  suggests  that  as  the  cyclone  becomes  more  intense  and  the  gcostrophic  winds  in- 
crease, downward  momentum  fluxes  begin  to  influence  frictional  forcing  more  than  they 
did  initially,  resulting  in  reduced  negative  frictional  forcing. 


2-1- li  surface  forcing  due  to  friction,  Case  I.  Forcing  in  K  100  k/tr]  dayx 
(dashed),  with  gcopotential  height  (solid)  in  m  and  observed  wind  vectors 
in  m  s~l.   Contour  interval  of  forcing  is  1. 


Fig.  23.        24-h  surface  forcing  due  to  diabatic  heating,  Case  I.    Contour  interval  is 
5. 


Forcing  due  to  diabatic  heating  at  the  surface  is  frontolytical  to  the  warm 
front,  with  the  values  of  forcing  larger  than  frictional  forcing.  Unlike  the  pattern  of 
frictional  forcing,  the  pattern  of  this  component  of  forcing  docs  not  resemble  the  surface 
heat  flux  pattern,  because  the  horizontal  gradient  of  the  heating  is  being  considered. 
What  is  evident  from  Fig.  23  h  that  the  surface  heating  gradient  tends  to  reduce  the 
low-level  temperature  gradient.  1  he  boundary  between  the  negative  forcing  region  along 
the  warm  front  and  the  positive  forcing  region  behind  the  front  implies  downward  mo- 
tion along  the  warm  front.  Consequently,  there  is  opposition  of  the  surface  warm 
frontogenesis  and  secondary  circulation.  The  diabatic  heating  component  dominates 
any  frontogenctical  effect  of  the  frictional  forcing,  giving  a  net  frontolytical  effect  of  the 
boundary  layer  processes  at  this  time. 


b.     Cross-Sections 

The  cross-sections  in  Figs.  24  and  25  depict  the  vertical  extent  of  these 
diabatic  forcing  terms.  Frictional  forcing  extends  up  to  900  nib  in  the  model,  at  which 
point  boundary  layer  stresses  arc  zero.  In  Fig.  24,  the  sense  of  flow  between  the  positive 
region  of  frictional  forcing  ahead  of  the  updraft  region  and  the  negative  region  beneath 
the  updraft  is  frontogenetical.  supporting  the  leading  edge  of  the  updraft  region. 
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Fig.  25.        As  in  Fig.  24,  except  for  diabatic  heating.    Contour  interval  is  5. 


The  boundary  between  the  negative  and  positive  regions  of  diabatic  forcing  beneath  the 
updraft  region  in  Fig.  25  is  frontolytical,  and  dominates  frictional  forcing  as  mentioned 
above.  However,  this  frontolytical  component  of  forcing  is  dominated  by  positive  values 
of  shear  and  confluent  geostrophic  dcformational  forcing  (not  shown)  in  the  boundary 
layer  ahead  of  the  front.  At  mid  and  upper  levels  diabatic  heating  generates 
frontogenesis  and  supports  forcing  of  the  updraft  region  throughout  its  vertical  extent. 
The  transport  of  heat  and  moisture  from  the  surface  to  the  mid  and  upper  levels  allows 
for  latent  heat  to  be  released  within  and  behind  the  warm  front,  which  enhances  the 
temperature  gradient  in  the  warm  air  sector.  Because  the  forcing  is  positive  it  also  en- 
hances the  secondary  circulation.  Geostrophic  deformation  is  also  positive  at  mid  levels, 
although  diabatic  heating  dominates,  especially  near  the  updraft  region.  The  diabatic 
forcing  is  most  susceptible  to  error  at  upper  levels  where  the  horizontal  temperature 
gradient  is  a  minimum,  so  a  further  discussion  of  upper-level  forcing  is  precluded,  and 
is  not  the  focus  of  this  study. 


3.     Case  2 

By  24  h  the  cyclone  in  this  case  has  deepened  to  9S4  mb.  Although  it  has  the 
strongest  warm  front,  it  has  about  the  same  accumulated  rainfall  as  Case  1  for  this  time, 
and  also  has  the  weakest  secondary  circulation  and  CCB.  The  patterns  of  confluent  and 
shear  geostrophic  deformation  are  similar  to  Case  1  and  are  presented  for  this  case  for 
reference. 

a.     Surface 

Figs.  26  and  27  depict  the  surface  forcing  due  to  confluent  and  shear 
geostrophic  deformation,  and  is  representative  of  the  typical  pattern  of  this  forcing  seen 
in  Case  1  after  24  h.  Fig.  2S  depicts  the  forcing  due  to  diabatic  heating  at  the  surface. 
With  the  absence  of  surface  heat  and  moisture  fluxes  in  this  case,  frictional  forcing  is 
negligible,  suggesting  that  surface  heat  fluxes  play  a  significant  role  in  driving  the 
downward  momentum  flux.  Although  there  are  no  surface  heat  and  moisture  fluxes  in 
this  ca^e.  the  pattern  of  forcing  due  to  diabatic  heating  at  the  surface  is  similar  to  that 
ol'  Case  1.  with  stronger  values.  This  suggests  that  either  the  errors  in  the  diabatic 
heating  calculation  are  relatively  large  or  that  there  are  significant  difference^  in  the 
boundary  layer  parameterization  in  the  model. 

The  sense  of  how  the  positive  and  negative  regions  of  forcing  enhance  or 
retard  the  secondary  circulation  i^  the  same  as  in  Case  1.  Again,  the  boundary  between 
the  band  o!  positive  forcing  near  the  v.  arm  front  in  Fig.  2S  and  the  broad  negative  region 
to  the  northeast  of  the  front  is  a  downward  motion  region.  Similarly  die  boundary  be- 
tween this  positive  region  and  the  negative  region  near  the  surface  low  is  upward.  As 
in  Case  i.  this  i^  Iron: 


Fig.  26.        24-h  confluent  geostrophic  forcing  at  the  surface.  Case  2.    Contour  inter- 
val is  5. 


Fig.  27.       As  in  Fig.  26.  except  for  forcing  due  to  shear.   Contour  interval  is  5. 


As  in  Fig.  26.  except  for  diabatic  heating.   Contour  interval  is  5. 


b.     Cross-Sections 

The  strength  of  the  vertical  circulation  is  weaker  in  this  case  than  in  Case 
1.  as  indicated  by  the  vertical  velocity  vectors  of  Iig^.  29  through  31.  Conllucnce  and 
shear,  in  Figs.  29  and  30  respectively,  are  strongest  in  the  boundary  layer  beneath  the 
updraft  region.  In  Fig.  31.  diabatic  heating  at  inid  levels,  in  the  form  of  latent  heat,  is 
not  as  strong  as  in  Case  1.  This  is  to  be  expected  since  the  only  direct  source  of  moisture 
in  the  model  is  from  the  PBL.  Not  only  is  there  no  frictional  confluence,  there  is  no 
moisture  transported  by  the  CCB  to  surface  low.  In  addition,  the  weaker  secondary 
circulation  in  this  case  supports  the  weaker  CCB  found  by  Nuss  (1989).  However,  this 
case  is  dominated  more  strongly  by  gcostrophic  deformation  at  the  surface  than  Case 
1.  which  supports  the  findings  in  previous  studies  that  gcostrophic  deformation  domi- 
nates frontogenesis  and  cyclogenesis,  since  this  case  has  the  strongest  warm  front  and 
eventually  the  deepest  cyclone. 
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Fig.  29.        24-h  cross-section  of  confluent  geostrophic  forcing.  Case  2.    Contour  in 
terval  is  5. 


Z4Z   XX  MARCH  IS 


SE2  EOT  Cvn 


Fig.  30.       As  in  Fig.  29,  except  for  forcing  due  to  shear.   Contour  interval  is  ; 
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Tig.  31.       As  in  Fig.  29,  except  for  diabatic  heating.    Contour  interval  is  5. 


4.     Case  3 

The  cyclone  in  this  case  has  developed  as  rapidly  as  the  Case  2  cyclone,  and  al- 
though the  surface  front  is  not  as  intense  as  in  Case  2  it  has  the  most  accumulated 
rainfall  at  24  h.  The  heat  (lux  distribution  at  24  h,  as  shown  in  Fig.  32,  is  characterized 
by  strong  upward  heat  fluxes  behind  the  cold  front,  positive  heat  fluxes  northeast  of  the 
surface  low  and  warm  from,  and  a  pocket  of  downward  heat  fluxes  near  the  trailing  edge 
of  the  center  of  the  warm  front,  similar  to  Case  1. 
a.     Surfat  c 

The  pattern  of  frictional  forcing  in  Case  3,  as  shown  in  Fig.  33,  is  quite 
different  to  that  of  Case  1,  in  that  almost  the  entire  region  cast  of  the  surface  low  is 
negative.  Within  this  negative  region  is  a  pocket  of  slightly  positive  frictional  forcing 
near  the  warm  front.  Again,  the  frictional  forcing  distribution  is  more  complex  than  the 
heat  flux  pattern,  indicating  the  significance  of  other  factors.  Whereas  in  Case  1  the 
region  of  positive  forcing  was  ahead  of  the  warm  front,  in  this  case  it  is  located  behind 


the  warm  front.  The  resultant  forcing  of  the  flow  between  the  positive  and  negative  re- 
gions does  not  seem  to  located  in  the  same  position  relative  to  the  warm  front  to  en- 
hance the  direct  circulation  cell  and  frontogenesis.  As  in  Case  1.  this  region  of  positive 
frictional  forcing  expands  in  si/e,  although  not  in  magnitude.  The  magnitude  of  negative 
frictional  forcing  does  increase  slightly,  although  it  decreases  in  area.  This  latter  result 
was  not  seen  in  Case  1,  and  suggests  that  since  the  two  cyclones  arc  of  equal  intensity 
at  this  time,  the  difference  is  due  to  the  surface  heat  flux  distribution. 

The  24-h  surface  diabatic  heating  in  this  case,  as  seen  in  Fig.  34,  is  similar 
to  the  other  cases,  but  with  a  more  narrow  band  of  negative  forcing  associated  with  the 
warm  front.  The  strength  of  the  gradient  from  the  negative  to  positive  regions  is 
strongest  in  this  case,  supporting  a  stronger  transverse  How  and  vertical  How  at  the 
boundaries.  Just  as  frictional  forcing  appears  to  be  less  frontogenctical  than  in  Case  1, 
the  locations  of  the  upward  and  downward  motion  regions  appear  to  be  less  front olytical 
than  in  the  other  two  cases.  An  analysis  of  the  cross  sections  yields  a  clearer  picture  of 
this  relationship. 


Fig.  32.        As  in  Fig.  21,  except  for  Case  3. 
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Tig.  33.       24-h  surface  forcing  due  to  friction.  Case  3.    Contour  interval  is  2. 


Fig.  34.        24-h  surface  forcing  due  to  diabatic  beating.  Case  3.    Contour  interval  is 


b.     Cross-Sections 

The  cross-section  of frictional  forcing  in  Fig.  35  reveals  that  the  region  of 
positive  frictional  forcing  is  locate  J  beneath  the  trailing  edge  (warm  sector  side)  of  the 
direct  circulation  cell  updraft  region,  with  negative  forcing  directly  beneath  and  ahead 
of  the  updraft  region.  The  resultant  vertical  motion  between  these  centers  of  forcing  acts 
to  oppose  the  vertical  component  of  the  direct  circulation  cell.  However,  because  the 
negative  frictional  forcing  ahead  of  the  front  is  strongest  in  this  case  compared  with  Case 
1,  one  can  imagine  a  stronger  transverse  (low  towards  the  front  and  cyclone  at  the  top 
of  the  PBL.  Thus  frictional  forcing  is  not  enhancing  the  updraft  region  but  is  enhancing 
the  transverse  component  of  the  direct  circulation  cell  at  low  levels. 

This  case  has  stronger  values  of  mid  level  latent  heat  than  the  other  two 
cases  at  this  time,  as  seen  in  Fig.  36,  and  indeed,  this  case  has  the  most  rainfall  accu- 
mulation at  24  h.    I  his  suggests  that  the  region  of  strong  upward  moisture  flux  northeast 


of  the  surface  front  has  influenced  the  latent  heat  of  the  mid  troposphere,  transported 
to  the  front  and  cyclone  by  the  influence  of  functional  forcins. 
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Tig.  35.       24-h  cross-section  of  frictional  forcing,  Case  3.   Contour  interval  is  1. 
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Fig.  36. 


As  in  Fig.  35,  except  lor  diabatic  heating.   Contour  interval  is  10. 


E.     TIME  SERIES  OF  FORCING 

The  results  at  24  h  arc  a  snapshot  at  one  time,  and  the  feedback  of  the  surface 
diabatic  effects  on  geostrophic  forcing  has  not  been  separated.  A  time  series  of  the 
contribution  to  total  forcing  by  each  component  at  low,  mid,  and  upper  levels  attempts 
to  separate  the  influence  of  each  forcing  term  as  they  evolve  throughout  the  simulation, 
forcing  was  measured  as  the  maximum  value  of  each  term  along  a  horizontal  line  across 
the  direct  circulation  cell  in  the  respective  cross-sections  of  each  case.  The  most  am- 
biguous term  at  low  levels  is  the  frictional  forcing  term.  Fractional  forcing  was  measured 
as  the  maximum  positive  value  near  the  leading  edge  of  the  updraft  region  of  the  direct 
circulation  cell,  since  this  would  enhance  the  updraft  region  in  this  location.  1  his  de- 
piction does  not  take  into  account  the  influence  of  the  negative  frictional  forcing  ahead 
of  the  warm  front.  Theoretically,  the  clockwise  circulation  around  the  negative  forcing 
center  is  acting  to  enhance  the  low-level  transverse  component  of  the  circulation  cell  at 
the  top  of  the  PBL  where  the  circulation  around  the  region  is  towards  the  front  and 


cyclone.  Forcing  due  to  diabatic  heating  was  measured  as  the  maximum  negative  value 
at  the  base  of  the  updraft  region,  since  it  was  generally  frontolytical  in  each  case  at  this 
level  and  location.  Confluence  and  shear  were  measured  as  the  maximum  positive  values 
in  the  direct  circulation  cell  within  and  ahead  of  the  front.  The  950  mb  level  was  chosen 
since  it  is  the  level  of  maximum  values  for  most  of  the  terms  at  each  time.  At  mid  levels, 
the  600  mb  level  was  chosen  to  represent  the  forcing  due  to  confluence,  shear,  and  latent 
heat  release.  Mid-level  forcing  was  found  to  be  a  maximum  at  this  level,  with  confluence 
and  shear  generally  centered  within  the  direct  circulation  cell,  and  latent  heating  located 
within  the  updraft  region.  Because  the  upper-level  diabatic  forcing  is  most  susceptible 
to  error,  the  maxima  were  located  at  300  mb.  where  the  center  of  the  upper-level  trans- 
verse component  is  located  in  each  case.  For  reference,  the  SLP  deepening  rate  is  also 
provided  with  the  950  mb  time  scries  in  each  case. 
1.     Boundary  Layer 

The  time  series  of  forcing  in  the  boundary  layer  for  Cases  1  and  2  is  depicted  in 
I  ig.  37.  and  for  Case  3  in  Fig.  38.  In  all  three  cases,  confluent  and  shear  geostrophic 
deformational  forcing  dominate  total  forcing  at  low  levels,  with  confluence  dominating 
shear  throughout,  the  simulations  except  for  Case  1,  where  shear  dominates  confluence. 
'I  otal  farcing  docs  not  peak  until  30  h  in  Case  3  and  36  h  in  Cases  1  and  2.  This  is 
consistent  with  the  decline  in  the  deepening  rate  of  the  cyclone  in  each  case  at  these 
times.  Surface  diabatic  heating  is  frontolytical.  initially  peaking  at  the  same  time  as 
geostrophic  forcing.  This  component  of  low  level  forcing  is  frontolytical  in  all  three 
cases,  and  reduces  total  forcing  by  about  20  "...  Unlike  geostrophic  deformational 
forcing,  diabatic  heating  at  the  surface  does  not  have  a  single  peak  and  period  of  decline. 
Frictional  forcing  is  not  a  significant  contributor  to  total  forcing  in  Cases  1  and  3.  but 
it  is  the  first  component  of  total  forcing  to  de\elop  in  Case  3.  In  all  three  cases. 
confluent  and  shear  geostrophic  deformational  forcing  and  diabatic  heating  do  not  de- 
velop at  low  levels  until  after  12  h  in  Case  3  and  IS  h  in  Cases  1  and  2.  The  earlier  de- 
velopment of  geostrophic  deformation  in  Case  3  is  presumably  due  to  the  influence  of 
the  frictional  forcing  that  was  present  from  the  beginning.  The  period  of  maximum  de- 
velopment of  geostrophic  and  diabatic  forcing  is  between  IS  h  and  2-4  h  in  all  three  cases. 
a;.  J  the  period  of  decline  is  between  30  h  and  36  h. 
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Fig.  37.        Time  series  of  forcing  in  the  PI3L.   (a)  Case  1,  and  (b)  Case  2.    Note  the 
scale  of  forcing  in  Case  2  is  twice  that  of  Case  1. 
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Fig.  38.        As  in  Fig.  37,  except  For  Case  3. 


2.     Mid  aiul  Upper  Le\els 

The  time  scries  of  total  forcing  at  mid  levels  is  depicted  in  figs.  39  and  4<>.  The 
time  scries  lor  upper  levels  in  depicted  in  figs.  41  and  42.  At  both  levels,  diabatic 
heating  dominates  total  forcing,  except  for  Case  3  where  geostrophic  deformational 

forcing  dominates.  Mid  level  diabatic  heating  in  the  form  of  latent  heat  release  is 
strongest  in  Case  1.  where  it  does  not  peak  and  level  off  until  42  h.  twelve  hours  after 
the  deepening  rate  has  begun  to  decline.  In  Cases  2  and  3  diabatic  heating  reaches  a 
steady  state  at  about  24  hours.  At  upper  levels,  there  is  some  contribution  to  total 
forcing  by  confluent  and  shear  geostrophic  deformation.  In  Cases  1  and  2,  diabatic 
heating  at  this  level  peaks  at  36  h.  then  declines,  whereas  in  Case  3  it  remains  steady 
after  24  h. 
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Fig.  41.       As  in  Fig.  37,  except  for  300  nib. 
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Fig.  42.       As  in  Fig.  38,  except  for  300  nib. 


F.     COMPARISON  OF  CASES 

1.  Boundary  Layer 

Case  2  is  similar  to  Case  1  during  the  first  18  hours,  alter  which  time  the  most 
significant  differences  are  the  rapid  rise  in  geostrophic  forcing  in  Case  2,  the  dominance 
of  shear  over  confluence  in  Case  1.  and  the  larger  surface  frontolytical  effects  of  diabatic 
heating  in  Case  2.  Comparing  Case  3  with  Case  1,  the  most  significant  differences  are 
the  stronger  frictional  forcing  in  the  first  12  hours  of  development,  the  rapid  rise  in 
geostrophic  forcing  similar  to  Case  2  lor  the  first  24  hour1;,  and  a  decline  in  the 
frontolytical  effects  of  diabatic  heating  after  30  hours. 

2.  Mid  and  I'pper  Levels 

The  most  significant  differences  between  Cases  1  and  2  arc  the  stead}-  rise  in 
diabatic  heating  in  Case  1  coupled  with  the  presence  of  geostrophic  forcing.  1  here  is 
no  significant  geostrophic  forcing  at  mid  levels  in  Case  2.  In  Case  3,  latent  heating  is 
evident,  but  it  not  the  dominant  mid  level  forcing  as  it  is  in  Case  1.  Both  confluent  and 
shear  deformational  forcing  dominate  diabatic  heating  in  Case  3.  At  upper  levels,  the 
only  significant  difference  is  the  absence  of  any  significant  geostrophic  deformational 
forcing  in  Case  2. 


V.     DISCUSSION 

A.  OVERVIEW 

The  basic  conceptual  model  of  frontogenesis  is  that  geostrophic  deformation  forces 
strong  frontogenesis  at  surface  and  upper  levels,  and  at  mid  levels  frontogenesis  is  usu- 
ally weak  and  forced  primarily  by  latent  heat  release.  Frontogenetical  forcing  can  be 
modified  by  small-scale  frictional  forcing  and  diabatic  heating,  with  their  relative  con- 
tributions to  frontogenesis  modified  by  surface  and  boundary  layer  processes. 

A  review  of  the  SLP  deepening  rates  from  Fig.  11  shows  that  for  the  first  12  hours 
each  case  deepens  at  the  same  rate,  indicating  that  the  surface  flux  distribution  peculiar 
to  each  case  doc;  not  have  an  apparent  effect  on  frontogenesis  until  12  h.  None  of  the 
cases  has  significant  geostrophic  deformation  enhancing  the  forcing  during  the  first  12 
hours,  suggesting  that  small-scale  processes  of  frictional  forcing  and  diabatic  heating  in 
the  boundary  layer  are  important  in  setting  the  stage  for  the  rate  of  development  and 
structure  of  the  cyclone  and  warm  front.  Once  cyclogenesis  has  begun,  forcing  due  to 
geostrophic  deformation  builds,  which  then  feeds  back  on  itself  and  dominates  the  forc- 
ing. In  later  stages,  friction  and  diabatic  heating  act  to  modify  the  geostrophic  defor- 
mational  forcing. 

B.  BOUNDARY  LAYER 

1.     Geostrophic  Deformational  Forcing 

In  their  two-dimensional  model  simulation  of  frontogenesis.  Baldwin  ct  al. 
(1984)  found  that  geostrophic  deformation  accounted  for  the  large-scale  features  of  the 
secondary  circulation  around  a  cold  front,  but  was  not  the  dominant  forcing  at  low  lev- 
els. Bond  and  Fleagle  (1985)  found  similar  results  in  their  case  study  of  an  actual 
oceanic  cold  front.  Baldwin  ct  al.  (1984)  found  that  shear  deformation  dominated  con- 
fluence in  their  model  simulation.  With  respect  to  the  forcing  of  a  warm  front,  this  study 
found  that  geostrophic  deformation  was  dominant  at  low  levels,  with  confluence  domi- 
nating shear  in  those  cases  with  the  most  intense  warm  fronts.  Mudrick  (1974)  found 
that  confluent  geostrophic  deformation  was  favorable  for  warm  frontogenesis.  with 
shear  geostrophic  deformation  favorable  for  cold  frontogenesis. 

As  mentioned  in  the  last  section,  confluence  dominates  shear  in  Cases  2  and  3. 
and  shear  geostrophic  deformational  forcing  dominates  confluence  at  all  levels  in  Case 
1.    Shear  deform  iti  >n  is  not  as  efficient  lor  frontogenesis  as  confluence  and  may  be  a 


partial  explanation  for  the  lower  intensity  of  the  warm  front  and  cyclone  in  Case  1.  The 
strong  frontolytical  effect  of  the  surface  fluxes  in  Case  1  to  reduce  the  temperature  gra- 
dient also  prevents  the  development  of  a  frontal  pressure  trough.  Lack  of  significant 
curvature  precludes  the  development  of  strong  confluence.  This  contrasts  with  Cases  2 
and  3.  winch  develop  stronger  temperature  gradients  and  more  of  a  pressure  trough 
along  the  front.  In  Case  3.  which  is  the  first  case  to  develop  significant  ge  atrophic  de- 
formational  forcing,  the  frontogenetical  effect  of  the  friction  apparenth  helps  to  con- 
centrate the  temperature  gradient  more  rapidly  in  this  case.  This  occurs  in  spite  of 
relatively  strong  frontolysis  due  to  surface  heating. 
2.     Frictional  Forcing 

Baldwin  ct  al.  (1984)  found  that  frictional  forcing  was  the  dominant  forcing  at 
low  levels,  and  Bond  and  Fleagle  (1^S5)  found  that  frictional  convergence  in  the 
boundary  layer  accounted  for  90%  of  the  observed  total  forcing  at  this  level.  The  results 
of  this  study  suggest  that  frictional  forcing  is  not  nearly  as  strong  a  component  of  forc- 

.  e  warm  front,  accounting  for  less  than  live  percent  of  forcing  at  low  levels.  In 
the  pattern  of  frictional  forcing  ahead  of  the  warm  front  is  not  as  consistent  in 
its  contribution  to  forcing  of  the  warm  front  as  it  is  to  the  forcing  of  the  cold  front  found 
in  these  other  studies.  It  i<.  not  known  whether  this  weaker  forcing  by  friction  in  the 
warm  front  studied  here  is  due  to  model  characteristics  or  some  fundamental  difference 
between  warm  and  cold  fronts.  The  stronger  frictional  forcing  and  more  intense  cold 
front  in  Case  3  compared  with  Case  1  suggest  that  frictional  forcing  can  greatly  enhance 
frontogenesis.     This  enhancement  is  partially  due  to  indirect  cHeus  on  other  more 

nt  processes. 
In  addition,  frictional  forcing  of  the  warm  front  is  not  consistently  predictable 
as  to  whether  it  is  frontogenetical.    As  shown  in  the  results,  sometimes  the  boundary 
between,  the  regions  of  positive  and  negative  forcing  was  such  as  to  oppose  the  forcing 
by  the  other  terms  and  sometimes  it  seems  to  enhance  the  forcing  by  the  other  terms. 
The  development  of  these  boundaries  seems  to  be  closely  related  to  the  position  of  the 
cycli  ne  and  warm  front  with  respect  to  the  surface  fluxes,  but  apparenth'  depends  on 
other  factors  as  well.     Surface  fluxes  remain  strongly  positive  behind  the  cold  front, 
whereas  the  flux  distribution  across  the  warm  front  is  confused  by  the  pocket  of  down- 
ward heat  lluxes  near  the  front  and  the  track  of  the  cyclone  relative  to  the  SST  gradient. 
There  is  the  suggestion  of  a  time  lag  between  frictional  forcing  and  its  effect  on 
gecurophic  deformational  forcing,  in  particular  confluent  forcing,  in  Case  3.   The  initial 
:    .        nent  of  confluent  forcing  is  preceded  by  12  hours  of  steady,  positive  frictional 


forcing  ahead  of  the  warm  front,  and  geostrophic  forcing  develops  6  hours  prior  to 
similar  forcing  in  the  other  two  cases,  which  have  little  or  no  frictional  forcing.  Simi- 
larly, the  decrease  in  the  rate  of  growth  of  the  confluent  (and  shear)  forcing  is  preceded 
by  12  hours  of  increasingly  negative  frictional  forcing  in  this  region.  The  reversal  in  the 
frictional  forcing  may  be  influenced  by  the  track  of  the  cyclone  relative  to  the  SST  gra- 
dient. During  the  first  24  hours,  the  cyclone  tracks  towards  the  region  of  strong  positive 
heat  fluxes  and  the  frictional  forcing  is  frontogenetical.  When  the  cyclone  begins  to 
move  into  colder  water  and  away  from  the  influence  of  the  positive  heat  fluxes,  the 
frictional  forcing  becomes  frontolytical. 
3.     Diabatic  Heating 

Baldwin  et  al.  (1984).  using  a  thermally  insulated  boundary  layer  in  their  model, 
found  that  boundary  layer  heat  lluxcs  had  a  minimal  effect  on  the  synoptic-scale  sec- 
ondary circulation.  Bond  and  Fleagle  ( 19S5)  included  boundary  layer  heat  lluxes  in  their 
study,  yet  found  similar  results  although  boundary  layer  turbulent  heat  lluxes  were 
found  to  be  frontolytical  in  the  boundary  layer.  Nuss  and  Anthes  (1987)  found  that 
boundary  layer  heat  fluxes  influenced  the  structure  only  below  850  mb.  suggesting  that 
surface  heat  flux  effects  are  confined  to  the  boundary  layer,  due  to  the  influence  of  the 
inversion  layer  above  the  PBL.  This  study  found  similar  results  in  that  forcing  due  to 
diabatic  healing  is  frontolytical  and  occurs  at  low  levels  only.  Its  magnitude  is  compa- 
rable to  that  of  shear  deformational  forcing  in  Case  2  and  confluent  forcing  in  Case  1. 
and  m  Case  3  it  is  about  half  that  of  shear.  Of  interest  is  the  presence  of  negative  forcing 
in  the  FBI  in  Case  2  despite  the  absence  of  surface  heat  and  moisture  fluxes.  This  sug- 
gests that  in  the  absence  of  differential  surface  fluxes,  there  is  surface  cooling  perceived 
by  the  model.  This  may  be  the  way  the  mode!  handles  convective  \ersus  non-convective 
precipitation,  or  it  may  be  a  result  of  the  computation  of  this  term  of  forcing  due  to 
advection. 

C.     MID  AND  UPPER  LEVELS 

1.     Geostrophic  Deformation 

As  with  studies  of  cold  front  forcing,  geostrophic  deformational  forcing  of  the 
warm  front  at  mid  levels  is  dominant  within  the  direct  circulation  cell  ahead  of  the  front, 
and  at  upper  levels  it  enhances  the  upper-level  transverse  component  of  the  secondary 
circulation.  However,  the  overall  values  are  typically  about  25"..  of  the  low-level 
geostrophic  forcing.  The  weakest  mid-  and  upper-level  geostrophic  forcing  is  in  Case  2. 
suggesting  that  mid-  and  upper-level  geostrophic  forcing  is  not  as  critical  to  the  devel- 


opment  of  the  cyclone  as  it  is  at  low  levels,  since  this  case  has  the  deepest  cyclone  and 
strongest  geostrophic  forcing  at  the  surface. 
2.     Diabatie  Heating 

Diabatic  heating  in  the  form  of  latent  heat  release  was  found  to  be  the  dominant 
forcing  mechanism  for  the  vertical  component  of  the  secondary  circulation,  with  positive 
centers  located  within  and  slightly  ahead  of  the  updraft  region  and  warm  front.  All  three 
cases  develop  some  mid-level  diabatic  forcing  due  to  latent  heating  in  nc  model.  The 
strongest  mid-level  latent  heating  developed  in  Case  1  and  it  continued  to  increase  until 
42  h.  This  case  was  the  weakest  cyclone,  which  indicates  that  latent  heating  and  a  strong 
vertical  circulation  in  the  warm  front  are  not  the  only  factors  for  strong  cyclogenesis. 
Case  3  developes  a  relatively  steady-state  of  latent  heating  by  24  h.  and  the  strongest 
vertical  motion,  however  the  diabatic  forcing  was  less  than  Case  1.  This  is  presumably 
a  feature  of  the  cumulus  parameterization  in  the  model.  Case  2  has  the  deepest  cyclone 
by  24  h  but  produced  the  weakest  mid-level  diabatic  forcing.  This  lack  of  latent  heating 
is  probably  a  reflection  of  the  lack  of  moisture  input  from  the  surface  in  this  case.  These 
results  suggest  that  details  of  the  cumulus  parameterization  and  moist  processes  in  the 
model  strongly  alter  the  intensity  of  the  mid-level  diabatic  forcing  of  frontogenesis  and 
the  secondary  circulation.  Since  the  diabatic  forcing  is  dominant  at  mid  and  upper  lev- 
els, the  circulation  at  these  levels  will  be  sensitive  to  model  errors  in  precipitation  proc- 
esses. 


VI.     CONCLUSIONS 

The  results  of  the  model  simulations  confirm  the  results  of  previous  observational 
and  numerical  studies  of  frontogenesis.  They  also  suggest  important  interactions  be- 
tween the  surface  heat  and  moisture  flux  distributions,  the  forcing  of  the  warm  front 
secondary  circulation,  and  cyclogenesis.  Geostrophic  deformational  forcing  was  domi- 
nant, but  evolved  differently  depending  upon  surface  friction  and  diabatic  heating  influ- 
ences. 

The  presence  of  frictional  forcing  enhances  the  low-level  transverse  component  of 
the  secondary  circulation,  which  in  turn  enhances  the  cold  conveyor  belt.  With  strong 
positive  heat  fluxes  in  a  favorable  location,  as  in  Case  3.  frictional  forcing  appears  to 
enhance  development  of  geostrophic  deformational  forcing  at  an  earlier  stage  of 
cyclogenesis.  which  also  enhances  the  secondary  circulation.  As  the  cyclone  and 
geostrophic  deformational  forcing  intensify,  the  surface  forcing  due  to  friction  acts  only 
to  modify  geostrophic  forcing.  The  development  of  positive  frictional  forcing  in  the 
latter  stages  of  cyclogenesis.  as  in  Case  1.  had  little  effect  on  the  deepening  rate,  except 
perhaps  to  maintain  the  geostrophic  forcing  at  such  a  level  as  to  yield  the  final  4  mb 
drop  before  cyclogenesis  finally  ceased. 

"1  he  effect  of  surface  heat  and  moisture  fluxes  and  their  distribution  on  forcing  of  the 
secondary  circulation  and  cyclogenesis  is  more  complex,  and  the  results  less  definitive 
than  with  frictional  forcing.  Surface  forcing  due  to  diabatic  heating  is  consistently 
frontolytical  in  each  case.  The  fact  that  it  is  a  strong  component  of  total  forcing  in  the 
absence  of  surface  heat  and  moistrure  fluxes,  as  in  Case  2.  suggests  that  either  the  error 
is  large  or  there  are  significant  boundary  layer  differences  in  the  model.  Early  surface 
relaxation  of  the  frontal  temperature  gradient  as  in  Case  1  appears  to  reduce  the  poten- 
tial deepening  rate  of  the  cyclone,  and  dominates  frictional  forcing. 


VII.     RECOMMENDATIONS 

The  results  suggest  areas  for  further  research,  including: 

Verily  the  magnitudes  of  the  various  forcing  terms  in  an  actual  case  study.  Because 
the  model  is  drier  than  real  oceanic  cyclones,  forcing  due  to  diabatic  heating  effects  may 
be  different. 

Verify  the  effect  of  surface  fluxes  and  boundary  layer  structure  on  the  frictional 
forcing  of  the  warm  front. 

Better  understand  the  role  of  various  model  parameterization s  involved  in  deter- 
mining convective  versus  non-convective  release  of  latent  heat.  This  may  help  to  explain 
the  presence  of  negative  surface  diabatic  forcing  in  the  absence  of  surface  heat  and 
moisture  fluxes. 
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